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Abstract 

When calculating higher terms of the e-expansion of massive Feynman diagrams, 
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are related to the derivatives of certain hypergeometric functions with respect to their 
parameters. We explore this connection and analytically calculate a number of such 
infinite sums, for an arbitrary value of the argument which corresponds to an arbitrary 
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1 Introduction 



In many cases, the results of analytical calculation of Feynman diagrams can be repre- 
sented as combinations of hypergeometric functions. However, the problem of constructing 
the e-expansion of hypergeometric functions within dimensional regularization [1] (where 
n = 4 — 2e is the space-time dimension) is not completely solved. Rather often, results 
for the terms of the e-expansion can be expressed in terms of polylogarithms [2] as well as 
Nielsen polylogarithms [3]. Recently it was demonstrated that new types of functions, har- 
monic polylogarithms [4] and multiple polylogarithms [5], appear in multiloop calculations. 
Furthermore, the nested sums [6] were proposed as a generalization of multiple polyloga- 
rithms, and the e-expansion for a large class of hypergeometric function was constructed [7] . 
In this paper we study the multiple inverse binomial sums defined as 1 

. . . 00 1 yj 

(i.i) 

where S a (j) = J2i=i^ a ^ s the harmonic sum 2 and u is an arbitrary argument. In what 
follows, we will also use z = \u as the argument of the occurring hypergeometric functions. 
For sums of the type (jl.ip . the weight J can be defined as J = c + J2k=i a kik + J2k=i bkjk, 
whereas the depth can be associated with the sum YX=i + J2l=i Jk- 

All multiple binomial sums can be presented in terms of function ip(z) = 4- lnT(z) 
and its derivatives by means of the following relation 

= {-l)\k - 1)! [C fc - S k (j - 1)] , k > 1, 

where ip( k >(z) is the fc-th derivative of the ^-function. In particular, for k = 1 we have 
= S\{j — 1) — 7s, where 7^ is Euler's constant. 
The sums (Jl.ip appear in the calculation of massive Feynman diagrams within several 
different approaches: for instance, as solutions of differential equations for Feynman am- 
plitudes [9] , through a naive e-expansion of hypergeometric functions within Mellin-Barnes 
technique [10], or in the framework of recently proposed algebraic approach [11]. Physi- 
cal applications include the one-, two- and three-loop massive Feynman diagrams with two 
massive cuts [12-17]. The case u — 1 corresponds to the single-scale propagator-type dia- 
grams [18]. Although there are many publications concerning harmonic series [19], only a 
limited number of results are available for the inverse binomial sums. Some particular results 
for u — 1, 2, 3 can be extracted from [16,20]. The sums with u — 1, 3 are expressible in terms 
of an "odd" basis [21], whereas the case u = 2 corresponds to an "even" basis [16, 22] 3 . 

1 Some particular results for the multiple binomial sums are presented in [8]. Those sums are defined 
similarly to Eq. but with ( 2 ?^j in the numerator, rather than in the denominator. 

2 Through the rest of this paper, the notations S a and S a will always mean S a (j — 1) and S a (2j — 1), 
respectively, even we do not mention this explicitly. When there are no sums of the type S a or Sb on the 
r.h.s. of Eq. we shall put a "— " sign instead of the indices (a,i) or (b,j) of X, respectively. 

3 The connection between "sixth root of unity" [23] and "odd" /"even" bases was discussed in [16,24]. 
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Several new results for a generalization of the sums (jl.lj) are presented in [25]. However, 
only for special type of sums the analytical results are available [20,26,27], 



where c > 2, 



6 = 2 arcsin (hVuj = 2 arcsin \fz , (1.3) 



2 sin — 
2 



Ls, (0) = Lsf } (5) (1.4) 



is the generalized log-sine function [2]. Some examples when such functions occur in the 
e-expansion of Feynman diagrams can be found in Refs. [21,28-30]. Usually, the occurring 
angles (jl.3j) . possess certain geometrical meaning [13]. 

The main aim of the present publication is the analytical calculation of inverse binomial 
sums (jl.lj) . including some relevant examples of physically important Feynman diagrams. 

The paper is organized as follows. In Section 2, employing the connection with the e- 
expansion of hypergeometric functions, we obtain analytical results for sums of the type (jl.lj) . 
valid for u < 4 (z < 1). Mainly the sums of the weights 3 and 4 are considered. In Section 3 
the analytical continuation to other values of u is constructed. Section 4 contains some 
applications of our results related to the e-expansion of Feynman diagrams, mainly two-loop 
master integrals. In Appendix A we briefly summarize the relevant properties of the har- 
monic polylogarithms of complex arguments and related functions. In Appendix B we show 
how certain identities between hypergeometric functions can be used to establish relations 
between the corresponding sums. In Appendix C we have collected explicit results for inverse 
binomial sums of lower weights. In Appendix D relations between binomial, harmonic and 
inverse binomial sums are explored, and analytical results for multiple binomial sums up to 
weight 3 are presented. Appendix E contains a realistic example of a physical application of 
the considered integrals. (Remember to change when appendices are ready!) 



2 Inverse binomial sums and hypergeometric functions 
2.1 The e-expansion of the 2 -Pi function 

Here, our analysis is based on comparing two representations of hypergeometric function 
whose parameters depend on e. One of them is the series representation in terms of the 
harmonic sums, whereas the second one is the exact result in terms of the functions related 
to the polylogarithms. 

In particular, we consider the 2 -Pi hypergeometric function of a special type, 



1 + a\e, 1 + (12S 
" + be 



2 



~ ^ {l + a l e) 3 {l + a 2 e) 3 {l + be) ] 
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where u = Az, (a)j = T(a + j)/T(a) is the Pochhammer symbol, and we have used the 
duplication formula (2/3) 2 j = 4?(/3)j(/3 + To perform the ^-expansion we use the well- 
known representation 



(1 + ae)j = j\ exp 



E 

L fc=i 



-as) 1 



k 



S k (j) 



(2.2) 



which yields 
2-^1 



1 + aie, 1 + a 2 e 
l + be 



2(1 + 2be) ^ 1 u j 



l + e 



u 



(A 1 + 6)51 - 265i 



+5 S 



26 2 (5 2 + 5 2 ) - 26(Ai + 6)5151 - \{A 2 + b 2 )S 2 + + 6) 2 5 2 



+ 6) 3 5 3 - 6(Ai + 6) 2 5 2 5i - |(Ai + 6)(A 2 + 6 2 )5i5 2 
+26 2 (A X + 6)51 (S 2 + 5 2 ) + b(A 2 + 6 2 )5 2 5i + |(26 3 - A 3 + 3AiA 2 )5 3 



-|6 3 (25 3 + 35 2 5! + 5 3 )] + G(e 4 )| , 

where A k = a\ + a 2 . Note that the coefficient of the e 3 term can be represented as 

l(Al + 6 3 )C + 2b 2 (A 1 + 6)Ci + 46 3 C 2 + \b{A^ - 6)(A X - 26)C 3 

+\{A\ - A 2 )(A ± + b) (5i5 2 - S 3 ) + ±6(A 2 - A 2 ) (s 3 - 25 2 5i) , 

where we have introduced the following combinations of harmonic sums: 

Cq = $i ~ 35i5 2 + 25 3 , 

C\ = S 3 — S\S 2 + Si(Si + S 2 ) — |5 2 5i + |5 2 5i , 

C 2 = jS\S 2 — IS* 3 + |5 2 5i — S 2 S\ — | (2S3 + 3SiS 2 + sfj , 

C 3 = S 3 — 2S 2 S\ — SiS 2 + 25 2 iSi . 



(2.3) 



(2.4) 



(2.5) 



To completely define the e 3 order of the expansion of 2 Fi, we need results for six combinations 
of sums, according to the number of independent combinations of the parameters A t and 6. 

Let us first consider the case 0<-u<4(0<-2<1). In this region the following 
parametrization can be used: u = 4 sin 2 | {z — sin 2 |), where < 9 < ir. In the rest of this 
paper, we will use the short-hand notation 

L e = In (2 cos f ) , l g = In (2 sin f) . (2.6) 

For a few special cases, the e-expansion of the 2 Fi functions is known. First of all, the 
following relation holds [31]: 



l+e,l-e 
2Fl { I 



sin 2 f = 



sin (e9) 
esin(6>) 



(2.7) 



Its expansion contains only the even powers of e. Using it, it is easy to get 



^ 1 V? 



$2 S^. 



|6> 3 tan | , 



^6> 5 tanf , 



(2.8) 
(2.9) 



etc. Then, there are some cases [16, 30] when an arbitrary term of the e-expansion can be 
calculated in terms of log-sine functions, 



F ( 1 ' 1 + £ 
2i H i + ie 



sin 2 f 



sin 2 f 



2 cos - ) 00 (Oc\i 

~ri E [ -T [Ls j+1 (7T - 0) - Ls J+1 (7T)] , (2.10) 



J'=0 



sin 



2 



sin #(2 sin fp^o J ! 
1 + e 00 e i 

-5:^(2^. 



(2 sin 



(2.11) 
(2.12) 



Moreover, for a more general case an integral representation can be obtained [16,31], 

2be— 2ae n 



1, 1 + ae 



l + be 



sin 



;i + 2fe) (2 cos §)' 
sin6 (2 sin I) 



/ a\2&£ / j\2ae— 2be 

d0(2sinf) (2 cos f J . (2.13) 



Any order of the e-expansion of Eq. (|2.13|) can be expressed in terms of the generalized 
"log-sine-cosine" function LsCjj(0) whose properties are listed in Appendix B of Ref. [16], 



Lscjj [9) 



In 



i-l 



2sinf 



ln J 



i-l 



2 cos I 



(2.14) 



Up to the level i + j = 5, only one independent function appears, Lsc 2)3 (6 l ), that cannot be 
expressed in terms of ordinary log-sine integrals. Using the definition ([2.14)1 . one can obtain 
the following representation: 

e 

Lsc 2 , 3 (#) = ^Ls 4 (20) - |Ls 4 (9) - f J d(f) In 3 (tanf ) . (2.15) 



The integral occurring in Eq. ()2.15j) can be expressed in terms of the inverse tangent integrals 
(see in Ref. [2]), 



1 r 



Tijv (z) = Im [U N (iz)] = — Li N (iz) - Li N (-iz 



2i 



Tijv (z) 



dx 



Tiiv-i (x) . (2.16) 
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For Ti^ (z), the following integral representation (see, e.g., in Ref. [28]) is useful: 



Ti/v (z) 



(N-iy. 



i + z 2 e 



(2.17) 



In this way, we arrive at the following representation of Lsc 2j 3(#) i n terms of the inverse 
tangent integrals (j2.16j) : 



Lsc 2>3 (0) = ^Ls 4 (26) - |Ls 4 (0) + 2Ti 4 (tan f) - 2 In (tan §) Ti 3 (tan § 
+ In 2 (tan f ) Ti 2 (tan f ) - \6 In 3 (tan f ) . 



(2.18) 



Note that Lsc 2j3 (7r — 6) can be reduced to Lsc 2j 3(0) by using Eq. (A. 27) of Ref. [16]. 

Using Eqs. ()2.10|) - (j2.13J) . we obtain the following results for the weight-2 and weight-3 
sums, in addition to Eq. (|2.8I) : 



00 1 v j 



2 tan 



00 1 u j 



00 1 u j 



tan 



00 1 7/J 



4 tan 



tan 



Ls 2 (-7T 
2LS 2 (7T 
LS 3 (7T 

| {5 [Ls 3 



0L fi 



(2.19) 
(2.20) 



) + Ls 2 (6) + Ok - 26L e 

) - Ls 3 (tt) - 2Ls 2 (tt-0) L e + + ^0 3 j ,(2.21) 

7) - Ls 3 (tt)] - Ls 3 (0) + |Ls 3 (26) - 2Ls 2 (0) L e 
+2Ls 2 (tt-0) Z fl - 8Ls 2 (tt-0) L e - 26l e L e + 40L 2 + ^0 3 }, (2.22) 



00 1 / 



tan 



I {6 [Ls 3 ( 



7T — 



Ls 3 (tt)] - 3Ls 3 (0) + 4Ls 2 (tt - 6) lg 



-8Ls 2 (tt - 6) L e + 2Ls 2 (0) l + Ls 3 (20) - 4Ls 2 (6) L 6 
+6l 2 9 - A6L e l e + 40L 2 + ^0 3 } , 



(2.23) 



where LS3 (7r) = — \itC,2- In fact, all Lsj (tt) are expressible in terms of the ^-function [2]. 

For the weight 4, the results for some combinations of sums involved in (J2.4j) can also be 
expressed in terms of log-sine functions, using Eqs. (|2.1(J|) - (j2.13j) . 



00 1 v j ( 
EtI7\— C ° = 8tanf{Ls 4 (7r-0)-Ls 4 (tt) 



1 v? 



-3 



Ls 3 (tt-0)-Ls 3 (tt)]l, + 3Ls 2 (tt-0) L 2 e -6L 3 e } 



(2.24) 



d = tan f {|Ls 4 (6) - §Ls 4 (26) - f [Ls 4 (tt - 6) - Ls 4 (tt)] 



+14 [Ls 3 (tt - 9) - Ls 3 (tt)] U - 2Ls 3 (9) k + Ls 3 (29) [L + l e ] - 14Ls 2 (tt - 9) L 2 9 
+2Ls 2 (vr -9)l 2 e - 2Ls 2 (29) L 9 l 9 - Ls 2 (29) L 2 - 29L e l 2 e - 29L 2 9 l 9 + A9L Z 9 ] , (2.25) 



00 1 v J < 

Em-^2 = tan §{±Ls 4 (29) - Ls 4 (9) + 4 [Ls 4 (tt - 9) - Ls 4 (tt)] - Ls 3 (29) (l 9 + L e ) 

-12 [Ls 3 (tt - 9) - Ls 3 (tt)] L e + 3Ls 3 (9) l 9 + 12Ls 2 (tt - 9) L 2 9 - 3Ls 2 (tt - 9) l 2 9 
-|Ls 2 (20) $ + Ls 2 (20) + 2Ls 2 (29) l e L e + 2^L + 29l e L 2 9 - \9l 3 - f9L 3 9 } ,(2.26) 

where Ls 4 (tt) = §7r£ 3 and the combinations of the harmonic sums Cj are defined in (J2.5|) . 
Finally, using f)2. 13j) the result for the sum involving C 3 can be expressed in terms of the 
Lsc-function f)2. 181) . 

E TiA = -4 tan §{2Lsc 2 , 3 (0) + 2 [Ls 4 (tt - 9) - Ls 4 (tt)] 
J =1 \ j ) 3 

+2 [Ls 3 (tt - 9) - Ls 3 (tt)] l e - 8 [Ls 3 (tt - 9) - Ls 3 (vr)] L 9 - Ls 3 (29) L 9 + 2Ls 3 (9) L e 
+8Ls 2 (tt - 9) L 2 6 - 4Ls 2 (tt - 9) L e l e + Ls 2 (29) L 2 9 - 29L 3 9 + 29L 2 9 l e } . (2.27) 

Some of the results ()2.19|) - (j2.27|) can be written in a slightly different form by means of 
relations [2] 

Ls 2 (9) = Cl 2 (9) , Cl 2n (tt + 9) = -Cl 2n (vr - 9) , Cl 2n+1 (vr + 9) = Cl 2n+1 (vr - 9) . 

However, in this way one cannot obtain results for the two remaining combinations 
in ()2.4|) . involving A 2 — A 2 = 2a^a 2 , since in Eq. f!2.13|) we always have a x a 2 = 0. An 
interesting relation between these two functions is obtained in Appendix B. It should also 
be noted that for higher hypergeometric functions other combinations of sums may arise at 
this level (see below). 

2.2 Expansion of higher functions 

Let us consider the hypergeometric function of the following type: 



p j | + b x e, . . . , | + 1 + a x e, . . . , 1 + a K e, 2 + d^ ... ,2 + d L e 

/>+l ; ' ' | + fie, . . . , | + fje, 1 + eie, . . . , 1 + e R e, 2 + c x e, ...,2 + c K+L - R - 2 e 



(2.28) 



where P = K + L + J — 2. Using the representation (j2.2j) its e-expansion can be reduced to 
inverse binomial sums (jl.lj) . The original hypergeometric function ()2.28|) can be written as 
(see details in Appendix B of Ref. [16]) 



P+l-rp 



{'l + ke}' 1 - 1 , {l + ai e} K , {2 + dis} L 
{§ + fre} J , {1 + e iS } R , {2 + Cl e} K+L - R ' 2 

2 Uf = \ L - R - 2 (l + c s e)U{ =1 (l + 2f k e) ^ 1 

" u nf =1 (l + d,e)n^(l + 2^) ft (f ) j^-^ 1 ' 1 ' 
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where u 
A = 



4z, 



exp 



E 

k=l 



-e 



(S fc T fc + 2 fc f45 fc + WfcT* 



1 - £ (W^r 1 + + 2^) + e 2 [\r 2 (W 2 + W 2 ) + W^- 1 (T X S X + 2^) 
+2T X U X S X S X + \T 2 S 2 + ±T 2 S 2 + 2 ([/ 2 £ 2 + f/ 2 S 2 



+ |j" 2 (W 2 + W 2 ) (T X S X + 2C/ 1 ^ ] 



-e 3 {|r 3 (2^3 + 3WW 2 + W{ 
+\W x] ~ l [t 2 S 2 + TfS\ + 4T X U X S X S X + 4 (c/ 2 s 2 + ulsi 
+\TlS\ + T 2 U X S X S X + \T X T 2 S X S 2 + T 2 U X S 2 S X + §T 3 S 3 
+2T 1 S' 1 (U 2 S 2 + E^S?) + I (U^Sf + 3i7i£7 2 ^ 2 + 2f/ 3 S 3 )} 



(2.30) 



and we introduced the constants 



Tk = Bk + Ck + Ek — Ak — Dk — Fk, 



N ^=E rf : 

Uk = Fk — Bk, 



W k = C k -D kl 



where the summations extend over all possible values of the parameters in Eq. (|2.28|) . We 
can see that for the general values of the parameters of the p+i-Fp function ()2.28|) we need all 
the occurring sums separately. However, for many applications, it is sufficient to consider the 
case J = 1 only. In this case, B k = and Uk = F k = f x . Substituting this into Eq. (|2.3U|) . we 
see that S 2 and S3 would only appear in combinations (S 2 + Sf) and (S'f + 3S X S 2 + 2S 
Furthermore, using the notations ()2.5|) . Eq. ()2.30|) in the case J = 1 can be presented as 



A 



j=i 



1 - s {W l3 - 1 + T X S X + 2/1&) + e 2 [lr 2 (W 2 + W 2 ) + W^ 1 (t x S x + 2f x S, 
+2f x T x S x S x + \T 2 S 2 + \TlSl + 2f\ (s 2 + S 2 



-r 2 



W 2 + Wf) [T x S x + 2f x S x 



-e 3 { Jr 3 (2W 3 + 3W x w 2 + wl 

+\W l3 ~ l [T 2 S 2 + T 2 S 2 + IfiT&S! + Afl (S 2 + Sf 
+\T X {t 2 + 3f x T x + 3/ 2 ) + \T X (T 2 - f x T x - f 2 ) S X S 2 + ±T 3 S 3 
+/i (T 2 + T 2 + 2f x T x + 2f 2 ) S 2 S X + 2 find - AffC 2 
-\h{T x + 2/0 (T x + 3/OCa} + 0{e A ) . 



(2.31) 



Eq. ()2.30|) also shows that we may need the sums with higher values of c (the power of 
which would come from 1/ j K ~ R ~ 1 in Eq. ()2.29|) and extra powers of 1/j in Eq. (|2.3U|) . 
Note that in the case K — R = 2, when we have the same number of and d{, we have 1/j 
factor in the sum in Eq. (|2.29J) . Moreover, when the parameters {2 + q} and {2 + di} are at 
all missing, we get no extra factors of 1/j in Eq. ()2.3()j) . for all Wk would vanish. 

For the sums of the type (jl.lj) with an arbitrary integer power c > 2 and < u < 4, the 
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following one-fold integral representation [27] is useful: 

" (c- 2)! sin * I ^ta^sm.Jj ^ ^ , (2.32) 

where /(j) stands for an arbitrary combination of the harmonic sums. In this manner, the 
original sum of the type with an arbitrary positive integer parameter c > 2 is reduced 
to a one-fold integral representation containing a sum with c = 1 in the integrand. Using the 
representation ()2.H2j) we can generalize our results ()2.8|) . ()2.9j) . (j2.19j) - ()2.27j) to the case of an 
arbitrary integer c > 1. However, only for special types of sums these results are expressible 
in terms of generalized log-sine functions (jl.4|) : 



00 u j 1 „ /- 2 (-2 



E ^TiA S * = -iE ? ,m 0™)^ Lag* (5) , (2.33) 



>'°(7) 



, : ■! \ ;• i , o ^( c 2 * 



00 u j 1 /„, „N 1 ^ 2 (-2 



where c > 2, < u < 4, and the angle 9 is defined in Eq. (jl.Hj) . 

Further results can be extracted from integral representation ()2.H2|) by using the integra- 
tion rules for Clausen's function [2] and the following relation: 

°r 1 
k / d0 ^Lsf (m0) = e k Ls { ; ] {mO) j-Ls^ (m6) . 



In this way, we obtain 



00 1 u j 

E/ITT-Si = 4C1 3 (tt - 0) - 2#C1 2 (tt - 0) + 3Cs, (2.35) 
00 1 7p - 

EyiA^i = -2Cl 3 (e)+4Cl 3 (7r-0)-2eCl 2 (7r-e)-eCl 2 (e) + 5C3,(2.36) 



Et4y^ 2 = ^ [Ls 3 (tt - 0) - Ls 3 (vr)] - 4 [Ls 2 (vr - 0)f + ^ 4 , (2.37) 



°° 1 7/ 3 ' 

E TiTYT^ 51 ^ 1 = 50[Ls 3 (7r-0)-Ls 3 (7r)]-0Ls 3 (0) + i0Ls3(20) 



-4 [Ls 2 (tt - 0)] 2 + - 2Ls 2 (tt - 6) Ls 2 (0) , (2.38) 
-4 [Ls 2 (tt - 6)f - [Ls 2 (9)] 2 + i0 4 - 4Ls 2 (tt - 9) Ls 2 (0) .(2.39) 
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However, other sums are not expressible in terms of the generalized log-sine functions. 
One needs to introduce a new function 4 , 



$(0) = j d(f) Ls 2 (0) In (2 cos f ) 



which obeys the following symmetry property: 

$(0) + $(tt - 0) = $(tt) + Ls 2 (tt - 0) Ls 2 (0) , 

where 

$(tt) = |ln 4 2-C 2 ln 2 2 + |C 3 ln2- fC 4 + 4Li 4 (i) = 0.64909... . 
The following sums of the weight 4 are expressible in terms of the function $(0): 



(2.40) 

(2.41) 
(2.42) 



Ls 



(i) 



7T 



Ls 



(i) 



7T 



It „(!) 



Ls 4 lj (20) + 2Ls 4 1} (0) - 40Ls 2 (tt - 0) Z e 



Cl 3 (vr)]/ 9 + 4$(vr 



Ls 4 i; (tt) - |Lsi 1} (20) + 4Ls 4 ij (0) + [Ls 2 



(i) 



00 1 u j 

-2vr [Ls 3 (tt - 0) - Ls 3 (tt)] + 8 [Cl 3 (tt 

-2tt [Ls 3 (tt - 0) - Ls 3 (tt)] + 8 [Cl 3 (tt - 0) - Cl 3 (tt)] l e - 4 [Cl 3 (0) - Ca] k 
-20Ls 2 (0) Z - 40Ls 2 (tt -9)l e + 4$(tt - 0) - 4$(tt) , 

where we have also used the following integral: 



4$(vr) , (2.43) 



(2.44) 



In (2 cos f ) In (2 sin f ) = - fLs^ (20) + ALb? 5 (0) + \ [Ls 4 1} (tt - 0) - Ls^ (tt 



-hr [Ls 3 (tt 



Ls 3 (tt)] 



Let us note that $(0) can be related to the real part of a certain harmonic polyloga- 
rithm [4] of complex argument, 



$(0) = <L0 2 (27r - 0) 2 - L Cl 3 (0) + C 3 m2 - H_ 1A0 ,i(l) + |[H_ 1AOil (e ie ) + H_ 1A0)1 (e 
where (for details, see Appendix A) 



(2.45) 



H-i,o,o,i0/) = [dx^P-, (2.46) 
J 1 + x 


H-i,o, ,i(l) = -^ln 4 2 + iC2ln 2 2-|C 3 ln2 + |C 4 -2Li 4 (|) = 0.33955... .(2.47) 

The function H_i ,0,0,1(2/) has a branch cut starting at the point y = — 1 (0 = ±7r), which is 
subtracted by logarithmic terms in (j2.45|) . so that $(0) is a smooth function of variable 0. 

Alternatively, instead of one may introduce the generalized Glashier function Gl4(#;l), see 

Eq. I|A.16(1 in Appendix A. 
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2.3 Further results for the sums 

For further investigation of the sums of the type (jl.ljl . a recursive approach appears to be 
useful. Some ideas applied below are described in the book [32]. Let us rewrite Eq. (jl.ljl in 
the following form: 

oo 

Kz£{fc± .(") = >-w;,(") = E " ! '/>:/,,.(./! , (2.48) 

where A = u\'"'.'^) and = f^'""^*J denote the collective sets of indices, whereas t/a-b-cU) 
is the coefficient of m j in Eq. (jl.ljl . i.e. the product of binomial sums divided by ( 2 ^) j c . 

The idea is to find a recurrence relation 5 , with respect to j, for the coefficients i]A;B;c{j) > 
and then transform it into a differential equation for the generating function S^ ; b ;c (m). In this 
way, the problem of summing the series would be reduced to solving a differential equation. 
Using the explicit form of T)A;B;c{j) given in Eq. (jl.ljl . the recurrence relation can be written 
in the following form: 

2(2j + 1)0" + l) c ~ W;cC7 + 1) - ./' '/.w,vt./) = r A Af) , (2-49) 
where the explicit form of the "remainder" ta-bU) is given by 

(?) r A ; B (j) = E[ [MJ - 1) +r afe f f[ K(2j - 1) + (2jr bi + (2j + l)~ 6 f 
fe=i i=i 

-f[[Sa k (j-l)rf[[S bl (2j-l)r . (2.50) 

fc=l z=l 

In other words, it contains all contributions generated by j~ ak , (2j)~ bl and (2j + l)~ bl which 
appear because of the shift of the index j. 

Multiplying both sides of Eq. (j2.49jl by u\ summing from 1 to infinity, and using the fact 
that any extra power of j corresponds to the derivative w(d/dw), we arrive at the following 
differential equation for the generating function T,a-b-,c{u): 

4 \/d\ c 2 / d \ c ~ * 

1 )( U T~) ^A;B;c{ U ) T, A -B;c{u) = ZVA-^-A 1 ) + R A,b{u) , (2.51) 

u / \ au J u \ au J 

where ?7^ ; b ;c (1) = \8 p0 and Ra-,b(u) = J2'jLi uj ' r A;B(j)- Using Eq. (j2.5(Jj) we obtain 

00 u j 1 

= ]^WjjZ> (2 - 52) 



00 w j 



1 



(2j) 61 (2j + 1) 



61 



(2.53) 



3 About application of the difference equations in the physical calculations we refer to [33] . 
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Up 



00 u> 



Sax _|_ ^2 _|_ 1 



Rai,a 2 ,a 3 ; — \U) 



S (?) 

h (?) 



h (?) 



j'a2 j'ai jai+a 2 



(2j) bl (2j + l)^ j 01 j ai (2j) 61 j'^(2j + l) 6 i 



(2.54) 
(2.55) 



<SW Sb 2 Sb ± 



+ 



5, 



6 2 



(2j) b2 (2j) fel (2j + l) fe2 (2j + l)^ (2j) 6l+62 
1 1 1 



jh( 2 j + l)b* ji»(2j + l)i* (2j + l) 



bi+b 2 



(2.56) 



U CL 2 1 



)'«2 



j'ai jd2+a:i jai+a-i, ja±+a 2 jai+a 2 +a 3 



In terms of the geometrical variable ()1.3j) . 



Furthermore, Eq. (|2.59|) can be represented as 



c-l 



tan • 



d6 



T, A . B . c (u e ) = tanf <t A ;b(0) 



where 



d 

dl 



(2.57) 



(2.58) 



.9/9 4 9/9 fl^- 

it = itg = 4 sin 5 , 1 = cot 7j , it— — = tan^ — , 

u 1 du z d9 

the differential equation (|2.51|) takes the following form: 

- si ^ 2 e (2 sin I cos 1^ - ij ( tan f ^) S A;jB . c (m ) = 5 p0 + R a ,b{u 6 ) . (2.59) 



(2.60) 



(7 

0"A;I?(#) = 5p0 + Ra,b(u 6 ) , r A; J B(6 l ) = S p0 8 + /# Ra-,b( U 4>) > ( 2 - 6 l) 



with = 4 sin 2 |. In particular, for c = 1 and c = 2 we have, respectively, 



SA;B;i(we) = tanf ct j 4; j b(6 i ) = tan- 



5 p0 # + dcj) Ra;b(u^) 



a a 

^A;B;2(ue) = J dcpa a -b{.(P) = Oa A ;B(d) - \5 p q6 2 - J cj)d(j) R a ,b(u 4 



Introducing lg = In (2 sin |) = | In ug, Eq. ()2.60|) can be rewritten as 



c—k 



\^f\ S A . B;c (M e ) = E A ;B;k(Ue) 



(2.62) 
(2.63) 



(2.64) 
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which is nothing but the differential form of the relation (|2.H2jl . The iterative solution of 
Eq. is 

Z A . B . c {u e ) = - 2^ — — l e ^A;B;c-i{u e ) + — — J d(j) — , (2.65) 

i=l l - Q ^ 

where = In ^2 sin f) = | lnw^. 

In particular, this solution allows us to formulate and prove the following important 
statement. If for some k the derivative T^A-.B-.c-kiue) is expressible only in terms of the 
powers of 6 and le, then the sum Y^a-b^uo) can be presented in terms of the generalized 
log-sine functions. Moreover, according to a statement proven in Appendix A.l of Ref. [16], 
the analytic continuation of any generalized log-sine function Ls^ (9) can be expressed in 
terms of Nielsen polylogarithms. 

Let us now explain how this general method works for specific sums of interest. Consider 

l - 00 1 u j } 

£ 3 ;_ ;1 (V) = H t^-v— = tanf / d</> -R 3; -(^) , 

i 1 I, J 3 o 
where Rs--(u) is defined in Eq. (J2.52|) and the result can be extracted from Eq. (jl.2j) : 

00 1 v j 

E TiTY- = 2 t Q 3 (*) + 9C ^ i°) ~ Cs] + , (2-66) 
Integrating over we obtain 

°° 1 7/-? r <i 

Em-* = tan f {6C1 4 (0) - e 2 C\ 2 (6) - 4#C1 3 (0) - 20Ca} • (2.67) 

For another sum, E 2 ' 1 C ;1 ( , u), we get 

1 1 _ 00 1 u j f 

S 2 ;ii- ; i(w) = 51 -p-y — ^1^2 = tan I y d</> i? 2 ,i;-K>) > 



where i2 2>1 ._ (u) is defined in Eq. (12341 . Using Eqs. d233), O and (l2~T)6l to calculate 
i?2,i ; -(u), and integrating over 0, we get 

00 1 77^ r 

= tanf{^Cl 2 (7r-^)-^Cl 2 (0)-40Cl3(7r-^)-40Cl3W 

-8CI4 (7T - 0) + 6CI4 (0) + 6(3 - \0 z L e } . (2.68) 
For the sum S 3 C 1 (u), this approach is also applicable. In this case we have 

s 3;-;i(«) = ~nj\~f s i = tan I / d< ^ > 

\ j ) 3 
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with #!,!,!;_(«) defined in Eq. (jZHZJ). Using Eqs. (j22U), and ijZTffijl to calculate 

and integrating over 0, we get 



00 1 v/- 7 ' r 

V^-Sj = tan|{6Cl 4 (0) - 24C1 4 (tt - - 129C\ 3 (n-9) - 40C1 3 (0) 
i ' (i ) ' 

-24L e [Ls 3 (?r - 0) - Ls 3 (tt)] + 8 [Ls 4 (tt - 0) - Ls 4 (tt)] - 2 C1 2 (0) 
+30 2 C1 2 (tt - 0) + 24C1 2 (tt - 0) L 2 e - 6 3 L g - 86L 3 e + 7( 3 6} . (2.69) 



Combining Eqs. (|2.67|) . (|2.68|) and ()2.69|) . we successfully reproduce Eq (|2.24j) . 
As a further example, let us consider the sum E^.^u). We obtain 

n 

00 1 v? r 

S 2;i;l(«) = 71T\~r S2Sl = taI1 f / d< ^ R 2;l{ U <t>) , 



3=1 



where i?2;i(w) corresponds to Eq. (|2.55j) . with some of the contributing sums having (2j + 1) 
in the denominator. Let us denote them as 



VA;B;c(j) ■ 



To calculate such sums, it is convenient to apply the differential operator that lowers the 
value of d, 



dw 



u 1/2 S^b^O) 



2 d 
cos f d0 



sin. T, A . B - C :d( u e) = T, A . B -c;d-i(ue) ■ 



[2.70) 



In particular, for d — 1 the sum on the r.h.s. is the standard sum (jl.lj) which is supposed to 
be known. Then the result for £a;-B;c;i(w) can be obtained by integration, 



^A;B;c;l(ug) 



Using the following decomposition 



2 sin 



d(fi cos f E y 4; jB;c (u ( , 



2 



a-l 



j«(2j + 



X / 8=0 



'6-1 + A 1 



6-1 



a — 1 + z 



(2j + 1) 



6-i ' 



we are able to rewrite T*A;B;c;d{u) as a linear combination of the sums, with one of the last 
two indices equal to zero: 



^; C » = E(-2)M d _[ S^oH + EH) 



i=0 



i=0 
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For the calculation of the T>A-,B;0;d(u) (with c = 0) we will apply d-times the procedure ()2.7()j) 

d / , \ 



2 d 



sin | £yi;B;0;d(we) = tan § — S A;B; i.o(«e) . 



d0 



This equation can be easily integrated for some particular cases 



OO 1 

E 



e 



+ i 



sin 



- 1 



sin 



Tii (tan f) - 1 



OO 1 

S(f)WTiT 
E 



— ^Ti 2 (tanf)-l 
sin h v y 



( 7 )(27+T) 51 = «co t | + — [Ls 2 (,-«)- 2 



oo j 

S(f)(2i + 1) 
^ 1 g 



ft(f)(2i + l) 



5, 



ft 



-^--20cot§-±0 2 + 4 
6 sin z z 



[2Ls 2 (vt-0)-20L, + Ls 2 (6) + 6l g ] 



(2.71) 

(2.72) 
(2.73) 
(2.74) 
(2.75) 



sin 9 



Si 



+ ^ cot f-— T Ti 2 (tanf)-l, 
sin 2 

-§0 2 + 4cot§ [Ls 2 (tt-0) -0L e ] -4 + 20 col £ + 



(2.76) 



3 



2 6 sin 9 



Ls 3 (7r-0)-Ls 3 (7r) + 0L 2 -2Ls 2 (vT-0)L e , (2.77) 



^ 1 v? 
fe( 2 /')(2j + l) 



S\S 2 



sin 

8-0 2 / e +20Cl 2 (tt-0)-20C1 2 (0)-4Cl 3 (vt-0)-2C1 3 (0)~Ca 



+ cot| 



±0 d -40+40L e -4Cl 2 (tt-6 

+-^{0 2 [Cl 2 (tt-0)-C1 2 (0)]-46/ [Cl 3 (tt-0) + Cl 3 (0)] 
sin (/ 

-8C1 4 (vt-0) + 6C1 4 (9) + 9( 3 -l9 3 L e } , (2.78) 

where Ti^r (z) is the inverse tangent integral (|2.1fij) . Together with Eqs. (|2.8jl and (|2.Hfij) . 
this provides us the result for R 2; i(u), 

R 2 ,i{u e ) = 4C1 3 (tt - 0) - Cl 3 (0) - 20C1 2 (tt - 0) +4( 3 + ±0% + 10 3 tan § + ^0 3 cot § . (2.79) 



Finally, using Eq. (J2.62)) and integrating over 0, we arrive at 

00 1 u j 



§ (?) ' 



5 2 ft = tan|{0 2 Cl 2 (vT-0) -40C1 3 (tt-0) -8Cl 4 (vr 
+|0 3 / 9 -i0% + 40C 3 } . 



Cl 4 (0) 



(2.80) 
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Using the results for the sums with c = 1 and applying Eq. (|2.32|) . we obtain the following 
results for the case c = 2: 



E 



1 ui 



S 3 = 9 2 C\ 3 (9)-69Ch(9)-12C\ 5 (9)-9 2 ( 3 + 12( 5 



[2.81] 



i =1 v J 
00 1 v j 



00 1 u j 



-\9 3 C\ 2 (tt - 0) + 29 2 C\ 3 (ir-9) + 9' 2 C\ 3 (9) + 80C1 4 (tt 



-Q9Ch (9) - 16C1 5 (tt - 9) - 12C1 5 (9) + ±6% - 3C 5 



(2.82) 



-|0 3 C1 2 (tt 



, , .., |0 3 C1 2 (9) + 29 2 C\ 3 (tt - 6 

+80C1 4 (vr - 0) + 0C1 4 (0) - 16C1 5 (tt - 9) + 2C1 5 (0) - 17Cs • (2.83) 



29 2 C 3 -\9 2 C\ 3 (9) 



We note that the results for other sums obtained in the previous sections can be also 
reproduced by this method. Moreover, it is possible to consider not only the combinations 
(S 2 + Sf) that correspond to the J = 1 case of Eq. (|2.3U|) (given in Eq. (|2.31j) ). but also 
the sums containing S 2 or Sf separately (which appear for J ^ 1). As an example, let us 
consider the sum 



00 1 u j 

S -;2;l(«) = E Tn7\ ~£>2 = tan- 



d(f> R^. 2 {u^) 



(2.84) 



where R^. i2 (u) corresponds to the case (|2.53|) . 
Substituting 

R-vM = zzn Ti 2 ( tan !) + 1^ 2 - 1 



sin 



into Eq. ()2.84j) and integrating over <ft, we arrive at 



1 7^YT^ = tan H 4Ti 3( tan I) + ^ 3 }- 



3 



(2.85) 



(2.86) 



Moreover, upon applying Eq. ()2.32j) we can perform another integration and obtain the result 
for c = 2, 

00 1 U j = . __. / „\ . r_. / „m2 



-S 2 = 49Ti 3 (tan f ) - 8 [Ti 2 (tan § 



+ 9^ 4 



^2.87) 



Using the connections between Ti 2 (tan |) , Ti 3 (tan |) and the log-sine functions (see Eq. (17) 
on p. 292 and Eq. (44) on p. 298 of Ref. [2]), we can get other representations, 



E-r^y— 5 2 = tan§{^ 3 + i#ln 2 (tanf) + 2 In (tan f) [Ls 2 (§) + Ls 2 (tt 



+2 



Ls 3 [it 



2J -Ls 3 (tt) 
16 



2Ls 



3 VI 



lLs 3 (9)}, 



(2- 



00 1 u j - 



Ls 3 [71 



20 



Ls 3 (1)1 + ±#Ls 3 (0) - 2 [Ls 2 (tt - §) + Ls 2 



(2.89) 



Note the appearance of log-sine functions of arguments | and 



3 Analytical continuation 

To obtain results valid in other regions of variable u (for u < and u > 4), we will construct 
the proper analytical continuation of the expressions presented in the previous section. For 
generalized log-sine integrals it is described in [16]. Let us introduce a new variable 



Me 



ln(— y — iaO) = Iny — i<77r, 



(3.1) 



where the choice of the sign o = ±1 is related to the causal "+i0" prescription for the 
propagators. For completeness, we also present the inverse relations, 



i - 



u — 4 



u- 



y d 
y- 



y 1 + \/^4 au 1+y dy 

and also expressions for (1 ± y) in terms of 8, 

1 - y = 2 sin § e'^^' 2 , 1 + y = 2 cos | e i,T9/2 . 



(3.2) 



(3.3) 



In terms of this variable y, the analytic continuation of all generalized log-sine integrals can 
be expressed in terms of Nielsen polylogarithms, whereas for the function <&(9) we get 

$(0) = C3ln2 + |C4-H_ 1A o,i(l) + |[H_i,o,o,i(l/) + H-i,o,o,i(y- i y 

- 1 l [u 4 (y) + Lu(y- 1 )]-^[H^ + y)-^y] [Li 3 (y) + Li 3 (r 1 )] > ( 3 - 4 ) 



with H_i )0 ,o,i(2/) defined in Eq. (I2~4T)|) . 

For the cases involving the inverse tangent integrals Tijv (tan|) and Tijv (tan|), the 
analytic continuation is straightforward (see Eq. ()2.16|) ). 



(tanf) 
Tijv (tan 



a 



— [Lijv (w) - Lijv (-w)] 



(7 



2i 



t [Li^ (o> 8 ) - Lijv (-^ s )] 



w = = — icrtan| . (3.5) 

1 + y 2 

u s = = — icrtan^ . (3.6) 

i + Vv 4 



Below we list the most complicated results, corresponding to the analytical continuation 
of the results obtained in Section 2: 



^ 1 v? 

he*)? 



Sl = -8S 1)2 (-y) lny + 4Li 3 {-y) Iny - 2Li 2 (-y) ln 2 y + 4 [Li 2 (-y)f 
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In 4 y + 4C 2 Li 2 (-y) + ( 2 In 2 y + 4( 3 In y + |C 4 , (3.7) 

E 71A -^^1 = -10Si >2 (-y) lny + S 1;2 (y 2 ) lny - 2S lj2 (y) lny + 4 [Li 2 (-y)] 2 
J=1 U J 3 

-2Li 2 (y) Li 2 (-y) + 3Li 3 (-y) lny - Li 3 (y) lny - §Li 2 (-y) In 2 y + ±Li 2 (y) In 2 y 

In 4 y + 6C 2 Li 2 (-y) - C 2 Li 2 (y) + |( 2 In 2 y + ^Cs lny + 5( 4 , (3.8) 

E TiA - + ^) = -12S 1)2 (-y) lny + 2S lj2 (y 2 ) lny - 6S 1)2 (y) lny + 4 [Li 2 (-y)f 
j=1 \ j) 3 

+ [Li 2 (y)] 2 - 4Li 2 (-y) Li 2 (y) + 2Li 3 (-y) lny - Li 3 (y) lny - Li 2 (-y) ln 2 y 

+ ±Li 2 (y) In 2 y + 8C 2 Li 2 (-y) - 4C 2 Li 2 (y) + ( 2 In 2 y + 8(3 In y + 10C 4 , (3.9) 

E 717Y- C i = YT^[ 3Li 4 (-?/) " 3Li 4 (1/) + 28S 1)3 (-y) - 2S lj3 (y 2 ) + 4S lj3 (y) 

-14S 2 , 2 (-y) + S 2 , 2 (y 2 ) - 2S 2 , 2 (y) + 28S 1)2 (-y) ln(l + y) - 2S lj2 (y 2 ) ln(l - y) 
-2S 1)2 (y 2 ) ln(l + y) + 4S 1)2 (y) ln(l - y) + 4Li 3 (-y) ln(l - y) + 2Li 3 (y) ln(l - y) 
-10Li 3 (-y) ln(l + y) + 4Li 3 (y) ln(l + y) - 2Li 2 (-y) ln 2 (l - y) 
-2Li 2 (y 2 ) ln(l + y) ln(l - y) + 12Li 2 (-y) ln 2 (l + y) - 2Li 2 (y) ln 2 (l + y) 

-21nyln 2 (l - y) ln(l + y) - 21nyln(l - y) ln 2 (l + y) + 4 lny ln 3 (l + y) 
+±ln 2 yln 2 (l - y) + 21n 2 yln(l - y) ln(l + y) - § In 2 y ln 2 (l + y) - § In 3 y ln(l - y) 
+| In 3 y ln(l + y) + f ( 4 - Cs ln(l - y) - 13& ln(l + y) + 7( 3 lny + f C2 In 2 y 
-C2 ln 2 (l - y) + 2( 2 ln(l - y) ln(l + y) + 8( 2 ln 2 (l + y) - 9( 2 lny ln(l + y)] , (3.10) 

°° 1 ?7 J : 1 — V r / \ / \ 

£ —y -C 2 = Y-^[2Si,3 (y 2 ) - 24S 1;3 (-y) - 6S 1>3 (y) + 12S 2 , 2 (-y) - S 2 , 2 (y 2 ) 

+3S 2 , 2 (y) - 2Li 4 (-y) + |Li 4 (y) - 24S 1)2 (-y) ln(l + y) - 6S lj2 (y) ln(l - y) 

-Li 3 (y) ln(l - y) + 2S lj2 (y 2 ) ln(l - y) + 2S 1)2 (y 2 ) ln(l + y) - 4Li 3 (-y) ln(l - y) 

+8Li 3 (-y) ln(l + y) - 4Li 3 (y) ln(l + y) + 2Li 2 (-y) ln 2 (l - y) - Li 2 (y) ln 2 (l - y) 

+4Li 2 (-y) ln(l - y) ln(l + y) + 4Li 2 (y) ln(l - y) ln(l + y) - 10Li 2 (-y) ln 2 (l + y) 

+2Li 2 (y) ln 2 (l + y) - § lny ln 3 (l - y) + 2 lny ln 2 (l - y) ln(l + y) 

+21nyln(l - y) ln 2 (l + y) - f lnyln 3 (l + y) - | In 2 y ln 2 (l - y) 

-2 In 2 y ln(l - y) ln(l + y) + 2 In 2 y ln 2 (l + y) + ^ In 3 y ln(l - y) - § In 3 y ln(l + y) 

In 4 y - |C 4 + 2C 3 ln(l - y) + llCa ln(l + y) - f Cs lny - 7( 2 ln 2 (l + y) 
+2( 2 ln 2 (l - y) - 2( 2 ln(l - y) ln(l + y) - f C2 In 2 y - C2 In y ln(l - y) 
+8C 2 lnyln(l + y)l , (3.11) 



oc 



E TiTY-^i = 4ff-i, 0> o,i(-y) + S 2 , 2 (y 2 - 4S 2 , 2 (y) - 4S 2 , 2 (-y) - 6Li 4 (-y) 
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-2Li 4 (y) + 4Si >2 (-?/) In y + 4S lj2 (y) In y - 2S lj2 (y 2 ) ln(y) + 4Li 3 (-y) ln(l - y) 
+2Li 3 (-y) In y + 2Li 3 (y) In y - Li 2 (y) In 2 y - 4Li 2 (-y) In y ln(l - y) 
-| In 3 y ln(l -t/) + ^ln 4 |/ + 2C 2 Li 2 (y) - |C 2 In 2 y + 2( 2 lny ln(l - y) 
+6C 3 ln(l-y)-3C 3 lny-4C4, (3.12) 

E TiA = ^-i,o,o,i(-y) + s 2, 2 (y 2 ) - 8S 2 , 2 (y) - 4S 2)2 (-y) - 6Li 4 (-y) 

• 2l.i, (y) - [Li 2 (y)] 2 + 4S 1)2 (-y) In y + 8S 1)2 (y) In y - 2S 1)2 (y 2 ) In y + ^ In 4 y 
+4Li 3 (-y) ln(l - y) - 4Li 3 (y) ln(l - y) + 2Li 3 (-y) In y - 4Li 2 (-y) In y ln(l - y) 
+2Li 2 (y) In y ln(l - y) - ±Li 2 (y) In 2 y - ± In 3 y ln(l - y) + 4( 2 In y ln(l - y) 
-Ca In 2 y + 10C 3 ln(l - y) - 5( 3 lny + 4C 2 Li 2 (y) - f ( 4 , (3.13) 
00 1 ?y J 1 — V r 

E TiTY-^i = h 48l M-J/) Ml + f ) - 48S i,3(-2/) + 24S 2 , 2 (-y) 
j =1 I -' 1 7 1 + v 



3 j 

-12C2 ln 2 (l + y) - 24 ln 2 (l + y)Li 2 (-y) + 24( 3 ln(l + y) + 24 ln(l + y)Li 3 (-y) 
-81nyln 3 (l + y) + 12( 2 lny ln(l + y) + 6 In 2 y ln 2 (l + y) - ln 3 yln(l + y) 
+^ In 4 y - |C 2 In 2 y + 3 In 2 yLi 2 (-y) + In 2 yLi 2 (y) - 5( 3 In y - 12 In yLi 3 (-y) 
-41nyLi 3 (y) + + 12Li 4 (-y) + 6Li 4 (y)l , (3.14) 



1 11^ 

E TiTT-^ = -12Li 6 (y) + 6 lnyLi 4 (y) - In 2 yLi 3 (y) - ^ In 5 y 

2 



j'=i 



+C 3 ln 2 y + 6C 4 lny + 12C 5 , (3.15) 

°° 1 M J 

E TiTT^i^ = -12Li 5 (y) - 16Li 5 (-y) + 61nyLi 4 (y) + 81nyLi 4 (-y) 
^' =1 I j J J 

- In 2 yLi 3 (y) - 2 In 2 yLi 3 (-y) + § In 3 yLi 2 (-y) + ^ In 5 y - i( 2 In 3 y 
-|C 3 ln 2 y-C4lny-3C 5 , (3.16) 
00 1 n? 

E TiA ~2 5,25,1 = 2Lis _ 16Lis (~^) _ m 2/ Li 4 (y) + 81nyLi 4 (-y) 
* =1 I i J J 

+| In 2 yLi 3 (y) - 2 In 2 yLi 3 (-y) + \ In 3 yLi 2 (-y) - \ In 3 yLi 2 (y) + ^ hi 5 y 

-|C 2 In 3 y - 2C 3 In 2 y - 8C 4 lny - 17C 5 , (3.17) 

where we have used the relations (|A.6|) and (|A.7|) . 

For the results involving Lsc 2j3 (9) we can use Eq. fl2.18j) to express them in terms of Ti 
function, and then employ Eq. (|3.5J) . For example, starting from Eq. (|C.5J) we obtain 

+2S li3 (y 2 ) - 8S 1)3 (y) - 48S 1)3 (-y) - S 2 , 2 (y 2 ) + 4S 2 , 2 (y) + 24S 2 , 2 (-y) 
+21n(l + y) [S li2 (y 2 ) - 24S lj2 (-y) + 10Li 3 (-y) - 2Li 3 (y) 
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81nyLi 3 (-y) 



+2 ln(l - y) [Si >2 (y 2 ) - 4S 1)2 (y) - 2Li 3 (-y) 

+2Li 2 (-y) [ln 2 (l - y) + 2 ln(l - y) ln(l + y) - 11 ln 2 (l + y) + In 2 y 
+2Li 2 (y) [ln 2 (l + y) - ln 2 (l - y) + 2 ln(l - y) ln(l + y)] + ± In 4 y 
+2 lny fln(l + y) ln 2 (l - y) + ln(l - y) ln 2 (l + y) - | ln 3 (l - y) - ^ ln 3 (l + y) 



-21n 2 yln(l - y) ln(l + y) + 5 In 2 y ln^(l + y) + |ln 3 yln(l - y) - fln 3 yln(l + y) 



2 i„3 



+C 2 31n 2 (l-y)-21n(l-y)ln(l+y)-131n 2 (l+y)-21nyln(l-y) + 141nyln(l+y)-f ln 2 y 



-HCs [3 ln(l - y) - 6 In y + 23 ln(l + y)] + f U 



(3.18) 



where u; = (1 — y)/(l + y), see Eq. f!3.5j) . We note that the analytic continuation of Lsc2,3(6*) 
can also be obtained directly, using the integral representation ()2.14|) . This procedure is 
described in Appendix A, Eqs. ()A.26j) - (IA.32l) . 

The analytic continuation of the sum involving £ 2 can be presented in terms of the 
variable uo s given in Eq. (|3.6|) . 

00 1 U 3 - r 1 r n2 

E727Y-^ = 21ny Li 3 (^)-Li 3 (-^) + 2 Li 2 (u 8 ) - Li 2 {-u s ) + ^hi 4 y. (3.19) 
' 1 [i) J 



The results for lower values of c can be deduced using 



d ~ v? 



For example, 



00 1 u 3 
00 1 U 3 - 



^ j c dw ^ J c+1 



4Li 3 (-y) - 2Li 2 (-y) In y - ± In 3 y + 3Cs + C2 In y 



-2Li 3 (y) + 4Li 3 (-y) - 2Li 2 (-y) lny + Li 2 (y) lny 



^ln 3 y + 2C 2 lny + 5C 3 



(3.20) 



(3.21) 



(3.22) 



The results for these two sums can be extracted from Ref. [14]. Further results for the sums 
are collected in Appendix C. 



4 Application to Feynman diagrams 

Below we present results for the e-expansion of one- and two-loop master integrals shown in 
Fig. 1. In the rest of this paper we use the notation u = p 2 /m 2 . 
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Figure 1: One- two- and three- loop diagrams considered in the paper. Bold and thin lines 
correspond to massive and massless propagators, respectively. 



4.1 One- loop vertex 



Let us consider a one-loop triangle diagram with mi = m 2 = m 3 = m, pf — p 2 — 0, with an 
arbitrary (off-shell) value of p\ = p 2 . Such diagrams occur, for example, in Higgs decay into 
two photons or two gluons via a massive quark loop. Following the notation of Ref. [10], 
we will denote this integral (with unit powers of propagators) as J 3 (l, 1, 1; m). According to 
Eq. (40) of Ref. [10], the result in an arbitrary space-time dimension n = 4 — 2e is 



J 3 (l, 1, 1; m) K , =0 = -\m 2 -%m 2 )-^V{l + e) 3 F 2 



1,1,1 + e 



5 2 

2' Z 



P 



Am 2 



(4.1) 



For this integral a number of one-fold integral representations are available, see Eqs. (3.10)- 
(3.11) in Ref. [16]. Expanding in e we get 



J 3 (l, 1, l;m) 



Ip2 =p 2 =0 



i7r 2 - £ r(l + e) 



[m 



2\-£ oo 



1 v? 



1 + eS x + \e 2 (S\ - S 2 



+±e 3 (Sf - 3SA + 2S 3 ) + ±e 4 (St - 6SlS 2 + SSiSs + 3S 2 2 - QS 4 ) +0(e 5 )] . (4.2) 
Substituting results for the occurring inverse binomial sums, we obtain 



J 3 (l, 1, l;m) 



=0 = - i7r 2 -r(l+e) 



[m 



2\-£ 
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+e 



AU 3 (-y)-2U 2 (-y) Iny - i In 3 y + ( 2 \ny + 3Csj 
2 [Li 2 (-y)] 2 + 2Li 3 {-v) \ny - 4S lj2 (-y) \ny - Li 2 (-y) In 2 ?/ 



+2C 2 Li 2 (-y) + |C 2 In 2 y + 2( 3 hi y - ^ In 4 y + |C 4 ] +0(e 3 ) } . (4.3) 

These results correspond to the analytic continuation of Eqs. (3.13), (3.14) and (3.16) from 
Ref. [16]. 

The e 3 term is Eq. (|4.2jl contains the same combination of sums as Co in Eq. (|2.5j) . 
However, the sum in Eq. ()4.2j) contains j 2 in the denominator, i.e., it corresponds to the 
case c = 2. Using the result ()2.24|) (with c = 1) together with ()2.32|) . we obtain a one-fold 
integral representation for the e 3 -term in Eq. (J4.2|) . Analyzing it, we see that we get only 
one new non-trivial integral, 



d0Ls 2 (0) In 2 



2sinf 



(4.4) 



while all other terms can be expressed in terms of known functions. For 9 = ^ the inte- 
gral (J4.4)) is connected with the new element %5 of the odd basis, whereas for 9 = | it is 
related to the new element x 5 of the even basis, see Section 3.3 in Ref. [16] and Eqs. (14)-(15) 
in Ref. [24]. 

At the same time, the e 4 -term of Eq. ()4.2|) (and all other even powers of e) can be 
calculated in terms of log-sine functions, so that their analytic continuation can be expressed 
in terms of Nielsen polylogarithms. Namely, for the £ 4 -term, combining Eqs. (3.14), (3.19) 
from Ref. [16] and Eq. (J2.34)) . we obtain the result for the sum 



£ (St - 6S?5 2 + 8^3 - 3S!) 

j=1 \ 3 ) 3 



with c = 2. Using Eq. ([3.20)1 . we also obtain the result for the case c = 1. 



4.2 Two-loop self-energy integral F 



10101 



This integral is a good illustration of the application of general expressions given in Section 3. 
The off-shell result for this integral in arbitrary dimension was presented in [34] (where it 
was called ^3, see Eq. (22) of [34]). For unit powers of propagators, the result reads 6 



2+ie, 



l-2e)F 10101 (p 2 ,m) 



;i- e *)(i+2 e ) 



4-^3 



l,l + e,l + e,l + 2e 



| + e,2 + e,2-£ 



Am 2 



1 



1.1+e.l+e 



2e(l + e 



■3^2 



2 • 



2 + e 



p- 



Am 2 



+ 



1 r 2 (i 



e 



2e r(l-2e) \ 
6 In given normalization each loop is divided by 7r 2_e r(l + e). 



m 2 Y 

2 3 ^ 2 

P J 



1,1,1 + e 
2 2 

2' Z 



Am 2 



.(4.5) 
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Expanding in e, we obtain 



p 2 (m 2 ) 2£ (l-2e)F wwl (p 2 ,m) 

r i 11 4 _ _ „ i 

0(e 2 



~ 1 u> [ 3 , . . 



-\ + — Si - - ln(- M )^ + \ ln 2 (- M ) - C 2 

J 2 J J 



= 6Li 3 (y) - 6Li 2 (y) In y - 2 In 2 y ln(l - y) - 6(3 

+e{28ff_ 1)0l0) i(-S/) + 7S 2 , 2 (V) - 28S 2 , 2 (-y) - 16S 2 , 2 (y) - 42Li 4 (-y) - 26Li 4 (y) 
+4 [Li 2 (y)] 2 + 16S li2 (y) lny - 14S 1)2 (y 2 ) lny + 28S lj2 (-y) lny + 18Li 3 (-y) lny 

+20Li 3 (y) In y + 20Li 3 (-y) ln(l - y) + 12Li 3 (y) ln(l - y) - 2Li 2 (-y) In 2 y - 9Li 2 (y) In 2 y 

-24Li 2 (-y) lny ln(l - y) - 4Li 2 (y) In y ln(l - y) - 2 In 3 y ln(l - y) + 6C 2 Li 2 (y) 

+4( 2 lny ln(l - y) - 12& lny + 24( 3 ln(l - y) - 9C 4 } + 0(£ 2 ) ■ (4.6) 

The result for the finite part coincides with [35] , whereas the result for the e-term is new. 
4.3 Two-loop sunset-type diagram Jon 

Let us consider sunset-type diagrams with two equal masses and one zero mass (see Fig. 1). 
The off-shell result for the sunset-type integral Jon with arbitrary powers of propagators has 
been obtained 7 in Refs. [34,36] by using the Mellin-Barnes technique [10]: 

T , 2 * , 2 , mi -, 2 r (^+"-") r (jK> r h±^zl) r (^+"-f 

J ii(cr, v u v 2 \P ,m) = [m ) -- - 



r(z/!)r(z/ 2 )r (f ) rfa + v 2 + 2a - n)r 2 (3 - § ) 

v F ( a;Z/H 
X 4-^3 n „ , 1 



cr, i^i + i/ 2 + cr - n, v 2 + a - |, v x + a - f 
2 , cr + |(z/i + z/ 2 - n), a + |(z/i + z/ 2 + 1 - n) 



^ . (4.7) 



4m 2 

For simplicity, in the definition of the integral Jon we will omit the arguments p 2 and m 

Jou(v, v\-,v-i) = Jon(cr, u u u 2 ;p 2 ,m). 

Let us remind that for the integrals Jon with different integer values of o and z/j there 
are two master integrals [38] of this type, Jon(l, 1, 1) and Jon(l, 1,2). However, two other 
independent combinations of the integrals of this type happen to be more suitable for con- 
structing the e-expansion, Jon(l, 2, 2) and [Jon(l, 2, 2) + 2J ii(2, 1, 2)] (see also in Ref. [39]). 
The latter combination corresponds to the integral Jon(l, 1, 1) in 2 — 2e dimensions [40]. In 
Ref. [18] Jon(l, 1,3) has been used as the second integral. To construct the e-expansion of 
the integrals Jon(l, 1, 1) and Jon(l, 1, 2) up to order e 2 , the integral Jon(l, 2, 2) should also 
be expanded up to e 2 , whereas [J n(l, 2, 2) + 2J 011 (2, 1, 2)] or Jon(l, 1, 3) up to the order e 
only. In Appendix C we give an example of a realistic calculation which demonstrates the 
required orders of the ^-expansion for the integrals involved. 

7 The imaginary part of the sunset diagrams in an arbitrary dimension is presented in [37] . 
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For Jon(l, 2, 2) we have obtained 

( 2\-l-2e 

Jon(l,2,2) = 



p ( l,l + e,l + 2e 
(l- £ )(l + 25) 3 2 1, § + e,2-e 



4m 2 



m 



2\-2e oo 



/r 



■E 



1 



+6° 



^ + i 











jl + e 


- + 55i - 


25i 






.J 






, 25 q2 
+ 


- 105151 




2 4- 25 2 + 25 2 



1,5c 2 c 1 25 o2 10 q q 5 q , 2 / q2 , q \ , 125 q3 ocrf o 



-f 5i5 2 + 4±5 3 + 55 2 5! + 105i 5 2 + 5 2 - f 5 3 + 35\5 2 + 25 3 



= 2- 



2\-2e 



(1 -y) £ y 



2e„.2eJ 1,„2 



ln^y 



+e [I In 3 y + C2 In y - 6 In yLi 2 (-y) - 4 In yLi 2 (y) + 3C 3 + 12Li 3 ( -y) + 6Li 3 (y)" 
+£ 2 [l2#_ 1A0 , 1 (-y) - 121nySi, 2 (-y) - 12 lnyS li2 (y 2 ) + 81nySi, 2 (y) - 12S 2 , 2 (-y) 
+3S 2 , 2 (y 2 ) - 12 ln(l - y)Li 3 (-y) - ^ hi 4 y - \( 2 In 2 y + 2( 3 In y + 6 In yLi 3 (-y) 



+8 lnyLi 3 (y) + 6C 2 Li 2 (-y) + 4C 2 Li 2 (y) + ±f- ( 4 + 18 [Li 2 (-y)] 
+12Li 2 (y) Li 2 (-y) + 4 [Li 2 (y)] 2 - f Li 4 (y 2 )] + 0(5 3 ) j , 

where the result for the e 2 -term is new. 
For another combination, we get 

,2\-l-2e 



(4.8) 



[m 



Joll (l,2,2) + 2/ on (2^,2) = -^^ 3 Fj 1 ' i 1 + + ; :i 1 +f 



2(m 2 )" 2£ g 



2 











55! - 25x" 


+ e 







Am 2 J 

f 5 2 - 105i5i - |5 2 + 25 2 + 25 2 



^fSi - 255i 2 5i - f 5i5 2 + 4±5 3 + 55 2 5i + 105i 5 2 + 5 2 



-| 5 1 3 + 35 1 5 2 + 25 3 



2(m2) ~ 2£ (1 y)1_2£ ^lny + 21n 2 y-C2-6Li 2 (-y)-2Li 2 (y) 



V 

+€ 



(i + y) 



l+6e 



24S li2 (-y) + 6S li2 (y 2 )-8S li2 (y) + |ln 3 y-4C 2 lny-llC 3 -24Li 3 (-y)-8Li 3 (y) 

24Li 4 (w) - 24Li 4 (-w) - 24Si, 3 (V ) + 12S 2 , 2 (y 2 ) + 64S lj3 (y) 
+96S 2 , 2 (-y) - 32S 2i2 (y) - 192S lj3 (-y) + 52Li 4 (y) + 48Li 4 (-y) 
-61ncu [2S lj2 (y 2 ) - 8S li2 (y) - 8S li2 (-y) + 7C 3 ] - 721nyLi 3 (-y) - 481nyLi 3 (y) 
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+6 In 2 cu [2Li 2 (y) - 2Li 2 (-y) - 3C 2 ] + 12 In 2 yLi 2 (y) + 18 In 2 yLi 2 (-y) 
+41nyln 3 cu + ^ln 4 y-5C 2 hi 2 y-26C3lny-57C4] +C>(e 3 )| , (4.9) 

where lu — (1 — y)/(l + y), see Eq. f!3.5|) . The e-term can be related to that of the result for 
J n(l, 1, 3) presented in Ref. [18], whereas the result for the e 2 term is new. 

4.4 Two-loop vertex diagrams 

Consider the two-loop vertex-type diagram P 126 given in Ref. [14] 8 (see Fig. 1), 




+\ \og 2 (-u) - S x log(-tt) - \S 2 - fS 2 + 4Si£i + 2y + .(4.10) 

It corresponds to one of the two-loop contributions to a boson decay into two massless 
particles, with a massive triangle subloop. Diagrams of such type have been intensively 
studied for Higgs boson production via gluon fusion [42]. 
Using our approach, we obtain 9 

(P 2 ) 2 ^i26 = ~ hi 2 y + ^ [2Li 2 (-y) my-4Li 3 (-y) + ln 2 y ln(l - y) -\ In 3 y-( 2 lny-3£ 
+BH- 1 ,o, , 1 (-y) + 2S 2 , 2 (y 2 ) - 8S 2)2 (-y) - 8S 2 , 2 (y) - 4Li 4 (y) - 12Li 4 (-y) 

-8Li 2 (y) Li 2 (-y) - 14 [Li 2 (-y)] 2 + 28S 1)2 (-y) lny - 10Li 3 (-y) lny 

+7Li 2 {-y) In 2 y - 4Li 2 (-y) Iny ln(l - - In 2 y ln 2 (l - y) - | In 4 y 

+ | In 3 y ln(l - y) + 2( 2 In y ln(l - y) - f C 2 In 2 y + 6(3 1b(1 - y) - HCa Iny 

-6Li 2 (-y) C2 - f C4 + 0(e) . (4.11) 

Here, the result for the finite part is new. Alternatively, results of such type can be obtained 
in a different way, using a technique based on Mellin-Barnes contour integrals [46]. 

Let us consider another two-loop vertex-type diagram shown in Fig. 1, K(vi, u 2 , a\, a 2 , 0-3), 
where two external momenta are on shell (p 2 = p 2 = m 2 ) and we also put p 2 = p 2 . Note 

8 The q 2 from Ref. [14] corresponds to our p 2 , whereas their z corresponds to our u. For P126 we keep the 
normalization used in in Ref. [14] for two-loop vertices: each loop integral is divided by 7r 2-£ m -2e e -TB£ , where 
je is Euler's constant. Note that it is different from the normalization used in Ref. [41], where each loop is 
divided by i-7r 2 ~ e r(l + e)^~ 2e , where \i is the scale parameter of dimensional regularization [1]. Furthermore, 
in numerical results presented in Table 1 of Ref. [41] an extra common factor m 4 r(l — 2e)/r 2 (l — e) was 
extracted. 

9 Recently this result was checked numerically by G. Passarino and S. Uccirati with the help of their 
approach [43]. After fixing a typo in an earlier version of our result, the results are in full agreement. 
The corrected result below was also confirmed (according to a private communication by R. Bonciani) by 
R. Bonciani, P. Mastrolia and E. Remiddi [44], using their approach [45]. 
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that K(l, 1, 1, 1, 0) appears as one of the master integrals in Ref. [45] (see Eqs. (102)-(106) 
of [45]). 

First of all, integration over the momentum of the massless loop (with powers of the 
propagators equal to 0\ and cr 2 ) can be easily performed, yielding the corresponding massless 
propagator to the power 0\ + o 2 — | , times a well-known factor containing T functions. The 
resulting integral is nothing but a one-loop triangle function, J 2 {y\, v 2) 01 + 02 + 03 — | ; m) 
in the notation of Ref. [10], with two external momenta (pi and p 2 ) on shell and z/ 3 = 
01+02+03 — 7*. Using the results of Ref. [10], one can see that in this on-shell limit the three- 
point integral J 2 reduces to a two-point function with two massive propagators (with 
the powers z/i and z/2), but with a shifted value of the space-time dimension (n — > n — 21/3), 



tt U3 r 2U3 



T(n -v x -v 2 - 2u 3 ) 
r(n -vx-v-i - ^3) 



J^\n-2v z -v x ,v 2 ) . (4.12) 



This is a generalization to arbitrary values of z/j of the relation found in Ref. [16, 47] for 
v x = v 2 = v z = 1. 



Using Eq. ()4.12|) . we arrive at the following result for the diagram K(ui, z/ 2 , 01, 02, 03; p 2 
(which we normalize by dividing each loop by i7r n//2 r(3 — §)): 

ir(z/i,z/ 2 ,0i,02,0 3 ;p 2 ,m) 

^4^14^3 r(z/!+Z/2 + 01 + 2 + 03-^)r (f-^l) T (f - ^) T (01 + 02~ 



[m 



2 \ i/i-H/ 2 +o-i+o-2+o"3^n 



r( ( x 1 )r( ( T2)r(z/ 1 +z/2)r(n- ( 7 1 -a 2 )r 2 (3 - f 



r (2n - v x - v 2 - 2o"i - 2o- 2 - 203) 
X tz ; — a-T< 



r ^-^1-^2-01-02-03 



Ui, v 2 , ui + u 2 + ai + a 2 + a 3 -n 
i, '- > 1 U"i + V2) 4(^ + ^ + 1) 




(4.13) 



where u = p 2 /m 2 . 

For v\ — v 2 — a\ — cr 2 = 1 and 03 = 0, the parameter z/ 3 in Eq. ()4.12j) gets equal to 
e, and the function corresponds to a one- loop two-point function in 4 — Ae dimensions (i.e., 
e -> 2e), 



#(1,1,1, 1,0;^,™) 



fm 2,-2 £ 1 r 2 (i- g )r(i + 2 £ )r(2-4e) 

1 1 2e 2 T(2-2e)T(2-3e)T(l+e) 2 1 




u 



4 • < 414 > 



As in the one-loop case (cf. Eqs. (C29)-(C30) of Ref. [47]), the 2 Fi function can be reduced 
to the type considered in section 2 by means of one of the Kummer relations, 




(l-4s) 2 F, 

All orders of the ^-expansion of the resulting 2 Fi function are given in (see also in Refs. [16, 
30]). For another special case, K(l, 1, 1, 1, 1), the parameter z/ 3 in Eq. ()4.12|) is equal to 1+e, 
and the function corresponds to a one-loop two-point function in 2 — Ae dimensions. In this 
case, we directly get the 2 Fi function of the type ()2.10|) . with e — > 2s. 
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Analytic continuation in terms of the Nielsen polylogarithms was discussed in section 2.2 
of Ref. [16]. In the case of K(l, 1, 1, 1, 0), we get 



#(1,1, 1,1,0;/, to) 



x 1 



TO 



2\-2s 



1 r 2 (l -e)r(l + 2e)r(l -4e) 



1 



' 2e 2 r(l-2e)r(l-3e)r(l + e) (l-2e)(l-3e) 

l-4e /i — 1_ 7/ U-4e oo 



x E EH 2 )"* S feji+1 _ fe _ p (-y) - (-lJ'S^+i-ifc-pr-y- 1 ) 
=0 z p. . 



fc=i 



(4.15) 



where the generalized polylogarithms Sfc iP (— y^ 1 ) can be expressed in terms of inverse argu- 
ment by means of the standard formulae given in Ref. [3]. In this way, we have obtained 
results for all coefficients of the e-expansion of #(1,1,1,1,0). The first three coefficients 
coincide with those given in Eqs. (104)-(106) of Ref. [45]. 



4.5 Three-loop vacuum diagram D4 

Consider a three- loop vacuum diagram with two different masses shown in Fig. (JTJ). Such 
an integral with equal masses (M = to) enters as a master integral in Avdeev's package [49] 
and MATAD [50]. Although there is no similar package for three-loop vacuum diagrams with 
two different masses, it is clear that such configuration (with unit powers of the propagators) 
will remain one of the master integrals in this more general case. 

This integral can be calculated by integrating over the momentum p the off-shell two- 
loop diagram iqoioi (also shown in Fig. with a massive propagator containing an arbitrary 
mass M (see Eq. (4.7) in [16]), 



D d (lAAAAA\u) 



17T 



d n p 



n/2 J p 2 _ M 2 



10101 



(p 2 ,m) 



(4.16) 



where u = M 2 /m 2 . The result of this integration can be presented in the following form: 
(m 2 ) 3£ (l- £ )(l-2 £ )£> 4 (l, 1, 1, 1, 1, 1; u) 

iF 3 



a 



1-6 



e(l-e)(l+e)(l+2e) 
Mr(l-g)r 2 (l+2e)r(l+3e) 
£(l+2 £ )(l+4e)r 2 (l+e)r(l+4e) 

l-e 



\A+eA+eA+2e 
|+e,2+e,2-e 



u 



2e 2 {l+e) 



3 r 2 




2e 2 (l + 2ef 



1 _ 2£ r(i- g )r(i+2 g ) 

U . ^77Z ; 3-^2 



Ae 2 T(l+e] 

^r(i-g)r 2 (i+2 £ )r(i+3£) 

'45 2 (l + 2 £ )(l+4e)r 2 (l+£)r(l+4£) 




l,l+2e,l+2e,l+3e 
3 
2 



| + 2£,2,2 + 2e 



3 r 2 



l,l+2e,l+2e 
l+e,2 + 2e 



1,1, 1+e 



+ 



r(i-£)r 2 (i+2e)r(i+3e) 

4e 4 r 2 (l+e)r(l+4£) 



l,l+2e, l+3e 
| + 2e,2 + 2e 



1 

ii 1 ' 



(4.17) 
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where, as before, u = Az = 4 sin §. All hypergeometric functions occurring in Eq. (|4.17|) 
belong to the type considered in section 2. Applying the results for the ^-expansion of these 
functions, we arrive at 

(m 2 ) 3e (l - e)(l - 2e)D A (l, 1, 1, 1, 1, 1; u) 

2Cs Q , , ^ «M 1 , 2 , 3 5 1 

= V - 9C4 + g(|)i-2? lnM + J - lnM -^- J - C2 

- ^ + - 1^ - + 1 + 5 a ) + 0(e)} 

= ^ + 2Lsi 1) (0)+8Z e [Cl 3 W-C3]-2^Ls 3 (0)-6 [Ls 2 (6>)] 2 + ^6> 4 - |C 2 ^ 2 - 9C 4 + 0(e).(4.18) 

Using analytic continuation described in section 3, we can also present this result in terms 
of the variable y defined in (13. 2|) , 

(m 2 ) 3£ (l - e)(l - 2e)D 4 (l, 1, 1, 1, 1, 1; u) 



2^3 1 2 2 

h 4 In it In y + In w 



6Li 3 (y) - 6 In yLi 2 (y) + | In 3 y - 3 In 2 y ln(l - y) - 6(3 



+4Li 4 (y) - 4S 2 , 2 (y) + 6 [Li 2 {y)f - 41n(l - y)Li 3 (y) + 121ny ln(l - y)Li 2 (y) 
-3 In 2 yU 2 (y) + 5 In 2 y ln 2 (l - y) - f In 3 y ln(l - y) + J In 4 y 

-12C 2 Li 2 (y) - 8C2 lny ln(l - y) + §C 2 hi 2 y + 4C 3 ln(l - y) + 3( 4 + 0(e). (4.19) 

One should remember however, that in this case the variable y becomes complex (although 
all imaginary parts should cancel for real masses), whereas Eq. ()4.18|) is explicitly real. 

As a non-trivial check on these results we consider two particular values, M 2 = m 2 and 
M 2 = 0. In the first case (9 = ~), we reproduce the known result (see in Ref. [23]) for the 
master integral D 4 = -D 4 (l, 1, 1, 1, 1, l;u)\ u ->i, 



,m 



2\3e 



4 - e)(l - 2e) D 4 (l, 1, 1, 1, 1, 1; u)\ u ^ = ^ - gC 4 - 6 [Ls 2 (f )]" + 0(e) . (4.20) 



In the second case, M 2 — (9 — 0, y — > 1), the result for the master integral Bm is 
reproduced (which was first calculated in Ref. [51]), 

(m 2 ) 3e (l-e)(l-2 £ ) D 4 (l, 1, 1, 1, 1, 1; «)U = (m 2 ) 3e (l-e)(l-2e)5 M = ^ - 9( 4 + 0(e) . 

(4.21) 

In a similar manner, analytical results can be deduced for other diagrams with two 
different mass scales, like D 3 and E 3 (for notations, see in Ref. [16]). To our knowledge, 
this is the first example of the calculation of the finite part of a three-loop vacuum diagram 
with six internal lines with two different mass scales. In previous publications [52], only the 
divergent parts of some three-loop vacuum integrals have been analyzed. 
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5 Conclusion 



In this paper, we have studied the multiple inverse binomial sums of the type ([l.ljh for 
arbitrary values of the variable u. Our analysis was based on the connection between these 
sums and terms of the expansion of certain hypergeometric functions with respect to the 
parameter e. Using known results for the hypergeometric functions, together with the inte- 
gral representation ()2.32j) . we have obtained a number of new analytical results up to the 
weight 4, Eqs. 02H2M22ZJ, (|2~55) - (j2~33jl and (I2~331) - (l2~m Moreover, in some cases like 
f!2.33|) - ()2.34|) the results can be obtained for an arbitrary weight. Constructing analytical 
continuation of the obtained results, we have expressed them in terms of the generalized 
polylogarithms ()3.7j) - (j3.13|) . In the cases considered, only one new function was needed, in 
addition to the basis of Nielsen polylogarithms. As such function one can take, e.g., the 
harmonic polylogarithm -ff-i,o,o,i( — y)> where the variable y is given in Eq. ()3.2|) . 

This approach allowed us to construct some terms of the e-expansion of the generalized 
hypergeometric function p + \Fp (I2.30|) and obtain new analytical results for higher terms 
of the e-expansion of some one- and two-loop propagator-type (see Eqs. ()4.6j) . ()4.8j) and 
(14. 9 j) ) and vertex-type (see Eqs. ()4.3j) and (|4.11j) ) diagrams depending on one dimensionless 
variable u = p 2 /m 2 . As a by-product, we have analytically proven the earlier published 
results for the three-loop vacuum integrals, the finite part of D4 and the e-part of D3 and 
E 3 , corresponding to the particular value of u — 1 (z — j), see Sections 4.3, 4.5 and 4.6 
in Ref. [16] and Refs. [21,23,61]. These integrals enter as master integrals in Avdeev's 
package [49] and MATAD [50]. Moreover, the developed technique is useful in the calculation 
of two-loop vertex-type diagrams [45,48]. 
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A Harmonic polylogarithms of complex arguments 

Here we collect some properties of the harmonic polylogarithms introduced in [4]. Let us 
define a w-dimensional vector a = {a,b), where a is the leftmost component of a, while b 
stands for the vector of the remaining (w — 1) components. The harmonic polylogarithms 
of weight w are then defined as follows: 




(Al) 
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where the three rational fractions f(a; x) are given by 

f(+l;x) = ^-, f(0;x) = -, f(-l;x) ' 



1 — X X 1 + X 

The derivatives can be written in a compact form 

^H 5 (y) = f(a;y)Hg(y). (A.2) 

Let us put y = e ld , so that 

6 

E s (e ie ) = H 3 (l) + i | d0 e 1 * f(o; e*) %(e^) , (A.3) 



where 



o 



e^+W/* 1 , 

^f(+i;^) = --^ = -- i-icot| y . 

2 sin | ^ v 
e^f(0; e i(f> ) = 1 , 

e^-lje*) = ? = - (l + itan f) . (A.4) 

2 cos | 2 v z ' 

Consider, for example, the harmonic polylogarithm #-i.o,o,i(y) given in Eq. (J2.46|) . Using 
the decomposition of Li 3 (e 1 ^ into the real and imaginary parts [2] and integrating by parts, 
we obtain expressions in terms of Clausen's and generalized log-sine functions, 



H-i A0 ,i(e ie ) = H_ lj0 , ,i(l) + ~ / # Li 3 (e^) (l + i 



S1I1 2 

cos I 



= H_ 1A0 ,i(l) - <^ 2 (2vr - Of + $(0) + In (2 cos §) Cl 3 (9) - Cs In 2 
+i In (2 cos f) Gl 3 (6) + ±i [Cl 4 (9) + Cl 4 (tt - 0)] 
-i [Gl 2 (9) Cl 2 (tt - 9) + !(tt - 0)C1 3 (tt - 0) - ±ttC1 3 (tt)] , (A.5) 

where 

Cl 3 (tt) = -|Cs , Gl 2 (9) =( 2 - \k9 + \9 2 , Gl 3 (9) = ±8 (tt - 0) (2vr - 0) , 

and $(0) is defined in Eq. fl2~HJ). 

In this paper, we have also used the following relations: 

#-1,0,0,1 = #-1,0,0,1 (-y) + Li 4 + Li 4 (y) - Li 3 (y) hiy - C2L12 (1 - y) 

+| In 2 2/Li 2 (y) + ± In 3 y ln(l - y) - §C 4 - iaTrf ( 3 ( A.6) 
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and 

H-i, QA i(y) + tf_ im i(-y) = Li 2 (y) Li 2 (-y) + ln(l + y)Li 3 (y) + ln(l - y)Li 3 (-y) . (A.7) 
The following two representations of -ff-i,o,o,i(l/) are a l so useful: 

H-wM = Li 4 (y) - |C 3 ln(l + y) + \ J ^ " , (A.8) 



00 



H-xmAv) = -|C3iu(i + y) + ^Ey 



(A.9) 



Instead of $(6*), one could also introduce another function. One of such possibilities is to 
consider a generalization of the Glashier function. Let us recall that the real part of Lij (e l9 ^J 
can be presented as 10 

Li^MH-tl^/^ W-^\n(l-2^ose + e) = i^ { ^ j . ° dd (A.10) 
lK J 2{j-2)\J £ V s s J \ Glj(6), 3 even v ; 

where Clj (8) and GL; (8) are Clausen and Glashier functions, respectively. In particular, 
GL, (8) is just a polynomial in 8. 

A possible non-trivial generalization of GL, (8) (for even j) could be 

Gl J .(g ;a ) = - — i / -^L i n i-2^ l n (l-2£cosfl + £ 2 ), j even, (A.ll) 



2(j-2)!7 £ + a 
so that Glj (8; 0) = GL,- (8). In particular 

Gl 2 (0;a) = -ln6 1n^+Li 2 (^,0)-Li 2 ^,0) , (A.12) 

where 



cosfl = -7= = = , (A.13) 



a + cos 8 
\/l + 2a cos # 
so that 

sin(# — 6 1 ) , > sin 6* . . 

a = — v ~ ; , b= Vl + 2acosd + a 2 = (A.14) 

sin sin 8 

Using Eqs. (17) and (18) on p. 293 of [2], we can see that the general result (jA.12|) simplifies 
in the case a — 1 (6 = |), 

Gl 2 (8; 1) = -\ Li 2 (cos 2 I) + |(tt - 8) 2 - \ In 2 2 . (A.15) 



3 We use the standard notation Lij (r, 0) = Re [Lij (re lS )] (see in Ref. [2]). 
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Then, let us consider 
Gl 4 (6 



(A.16) 



For 9 = we get 

Gl 4 (0; 1) = C/3,1 - C4 , (A.17) 
where L^i is the alternating two-fold Euler sum considered in [22,23] (see also Ref. [53]), 

U 3A = -2Li 4 (±) + \U~T2 ln * 2 + IC2 In 2 2 . 



For general one can see that GI4 (9; 1) is related to the real part of a harmonic poly- 
logarithm, 

Re#_ 1A0 ,i(e ie ) = GI4 (9) - GI4 (0; 1) - |Ca In (2 cos f ) . (A.18) 
In particular, the function $(6 I ) can be presented as 

*(0) = -GI4 {6; 1) + Gl 4 (0; 1) + |Csln2 - ±9 2 (2n - 9) 



where Cl 3 (tt) = — |Cs- 

Using the symmetry property of $>(9), Eq. (j2.41|) . we get 



[Cl 3 (O)-Ch (7r)]ln (2cos - 2y . 

(A.19) 



Gl 4 (0; 1) + GI4 (tt - 9; 1) 



11 



C4 - jz0 2 ( 



■18 



TT 



Cl 2 (0) Cl 2 (tt 



-[Cl 3 W-CI3 (7r)]ln (2 cos I) 
-[Cl 3 (tt - 9) - Cl 3 (tt)] In (2 sin § 



For special values of 9, Eq. (jA.20|) yields 

GUM) = 



GU ( f ; 1 



-C/ 3l l + fC4-K3 ln2 > 

HC4 - §Csln2 - G 2 , 



Gl 4 (f;l)+GU(f;l) = %U - iCs In 3 - f [d 2 (f ) 



(A.20) 

(A.21) 
(A.22) 

(A.23) 



where G? is the Catalan constant. Moreover, with the help of PSLQ algorithm [54] one can 
obtain results for GU ( |; lj and GU (^-; lj separately, 



GU(f;l) - ffC4 - MCsinS - I 

Gl 4 (f;l) = 



Cl 2 f5 



216^ 4 



Cl2 f 



+ I^ Ls 3 (f ) - §Ls?> (f ) , (A.24) 
17TLS3 (f ) + fLsf J (?) . (A.25) 



We would also like to discuss the analytical continuation of Lsc 2i3 (#) (see Eq. (12.14)1 ). 



Lsc 2 , 3 (0) 



ln 



2 sin I 



In 1 



2 cos ■ 
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If we introduce a variable z = e lcr< ^, we see that 

In (2 sin f) <->ln(l-3) - §lnz+±i<77r, In (2cos|) ln(l + z) - \\n.z, (A.26) 
and the analytical continuation of Lsc 2j3 (#) is given by 



d2 



iaLsc 2i3 (#)= / — {ln(l-z)-iln2 + |iavr}[ln 2 (l+z)-ln2ln(l + z) + iln 2 



with y defined in (|3.1|) . This integral can be calculated by using the following relations 



/dz 
— \n(l + z) ln(l-z) In z 



J^\n\lTz) ln(l±z) 



S 2 , 2 (y) + S 2 , 2 (-y) - iS 2 , 2 (y 2 ) 

+ lny[iSi j2 (?/ 2 ) - Si >2 (2/) - Si, 2 (— 2/) 

21n(l±y)S 1 , 2 (±y)-2H_ 1A i,i(±y), 



where 



H_i, ,i,i(y) = / 7-^- Si >2 («) 
J 1+2; 


is another harmonic polylogarithm [19]. We note that 

y , 

H_ li o, 1 , 1 (±y) = iS li3 (y 2 )-S 1 ,3(^)+ln(l±y)S 1 , 2 (± 2 /)±^ | In 3 (- 



1-2 



+ -2 



H_i, ,i,i(y) + H_i A i,i(-2/) = iSi i3 (y 2 J - Si, 3 (y) - Si, 3 (-y) 

+ ln(l + y)S 1 , 2 (y)+ln(l-y)S 1 , 2 (-y) . 

Therefore, the analytic continuation of Lsc 2j3 (6 l ) reads 

icxLsc 2 , 3 (#) = S 2 , 2 (y) - |S 2 , 2 (y 2 ) - 2S 1)3 (y) + ±S lj3 (y 2 ) + ±Li 4 (y) + Li 4 (-y) 
+ In y [§S li2 (y 2 ) - S 1)2 (y) - |Li 3 (y) - Li 3 (-y)' 

"dz 



(A.27) 
(A.28) 

(A.29) 
(A.30) 
(A.31) 



+ In 2 y ±Li 2 (y) + ±Li 2 (-y) + i In 4 y - ^4 - | / — ln : 



1 - 2 

1 + 2 
U 

-iavr {S 1|2 (-y) - |Li 3 (-y) + §lnyLi 2 (-y) + ^ In 3 y - |C 3 } • (A.32) 

Note that the integral occurring in Eqs. (|A.30|) and (|A.32|) is directly related to the 
integral in Eq. ([2.150 . and it can be calculated in terms of polylogarithms, 

y 



dz 



In 3 



1 - z 

1 + 2 



-f C 4 + 6 [Li 4 (u) - Li 4 (-co)) -Q\mu [Li 3 (u) - Li 3 (-u)} 

+3 In 2 u [Li 2 (u) - Li 2 (-w)] + In 3 u In y , ( A.33) 
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where to = (l — y)/(l+y) (see Eq. (I3.5|l ). Therefore, H_i j0) i ) i(±|/) reduces to polylogarithms, 
as well as the analytic continuation of Lsc 2j 3(#). We note that the imaginary part on the 
r.h.s. of Eq. (|A.32|) gets cancelled in the results for the corresponding sums. 
The following relations are useful for the transformations: 

Li 2 (u) - Li 2 (-w) = -Li 2 (y) + Li 2 (-y)- In In y+§C2 , (A.34) 
Li 3 (u) - Li 3 (-w) = ^1,2 (y 2 ) - 2S 1)2 (y) - 2S 1>2 (-3/) - In a; [Li 2 (y) - Li 2 (-</)] 

-iln^lny + fCaln^ + ICs • (A.35) 



B Additional identities between inverse binomial sums 



As we have seen, the e-expansion of the hypergeometric functions produces series of the type 
(II. lj) . One can use certain properties of hypergeometric functions to get relations between 
the sums (jl.lj) involved in the ^-expansion. 

Let us consider Eq. (B.18) of Ref. [16] (which follows from Eq. (22) on p. 498 of [31]), 



3-^2 



1 + aie, 1 + a 2 e, 1 + \{a x + a 2 )e 
I + |(ai + a 2 )e, 2 + (a x + a 2 )£ 



2^1 



1 + ~ai£, 1 + \a 2 e 
| + |(ai + a 2 )e 



(B.l) 



It reduces the given 3 F 2 function to a square of the 2 Fi function. Substituting the e- 
expansions of the ^F 2 and 2 F\ functions into (jB.ljl . we obtain the following relations between 
the sums: 

„ — , s 4 — It r _._ , 



[U 



2u 



u 



SlL ;2 («) - S_;J. 2 (u) - S_:_. 3 (m) 
, s 4- w 



4 - w 



J 2;-;2 



4u 



s_L.i(M)si;_ ;1 (w) 



(B.2) 
(B.3) 
(B.4) 



>:?; ;2 («) - 2Eg. 2 («) + e:;?. 2 («) + s:g. a («) - 2E};:. 3 («) + 2s:;} ;3 (n) + 2s: 
1 " x(«) - 2Eg ;1 («) + e:; 2 ;1 («) + >: ; 2 » - |s 2 ;: ;1 



4 (M 



(B.5) 



with m = 4z, so that 



4 - u 



2 ' 



li 2! 

Another interesting relation, Eq. (B.19) of Ref. [16] (which follows from Eq. (20) on p. 498 
of [31]), reads 



3^2 



1 + axe, 1 + a 2 e, 1 + \{a x + a 2 )e 
\(<xi + a 2 )e, 1 + (a x + a 2 )£ 



2 1 2 

xfl 



/ 1 + |ax£, 1 + |a 2 £ 
2 1 ^ | + \{a 1 + a 2 )e 



a\(x 2 8 z 



2 [1 + (oi + a 2 )e] 



3^2 



1 + |oi£, 1 + |a 2 £, 1 
f + |(ai + a 2 )e, 2 



(B.6) 
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In this way, we get the following relations: 



§E_:_ ;1 (?/)E_:_. 2 U' s 

2' 

:2t " ! - 



^2;-;l( U ) — 3^— ;— ;1 V"-^^— ;— ;2 V^^ ) 



.(B.7) 



We have checked that all above equations are satisfied by the explicit results for the sums 
listed in this paper. 

Let us also present some relations for higher-order sums: 



S 3!-;2W + E 2 L. 3 («) ~ ^H-A U ) + Sgl} 2 («) 



A-u 
Au 



{ S 2i-;lH 



+ [l^lM - S 2;i;-;l( M ) + | S 3!-;l( M )]} 



S 2i-;2(^) - E^L.^) 



a(«) 



J 4|-;l 



A-u 
16m 



E 2i-;lN 



El:: ;1 (u) 



S 2;-;l( M ) - S 4;-;l( M ) 



+ 3E 2 :_ Ju) 



(B.8) 



(B.9) 









2 - 


2 


Eii: 







+6E 2 ;I ;2 ( M ) 

.2:1 



S 5|-;l( W ) ~ + 4S -i-;2(«) S -i-;lH [ S li-;2(«) - 



+45E^ ;1 ( M ) - 45Eji;:. 1 («) + 45E^::.» - 45£# ;1 (u) + 93^; 2 ;: ;1 (u) - 93E£: ;1 (u) 



-d,l;- 



1:1 



.1,1;- 



+4s^;: ;1 ( u )e: 



(u) - 3E 



3;-;2 



(«)e: 







(B.10) 



where Eqs. ffB~8j) and (IB~9j) follow from Eq. (IB~T| . while Eqs. (lB~9j) and (lB~T0|) follow from 
Eq. (|B~6|) . 

One more relation can be derived from Eq. (47) on p. 456 of [31], 



1 + a x e, 1 + a 2 e 

3 + 1 
2^2 



3 + §(% + a 2 )e 



sm 



COS6 1 2 iq 



1 + \a\e, 1 + \a 2 e 
§ + |(cti + a 2 )e 



sin 2 . 



(B.ll) 



We can compare the e-expansions of these 2 F\ functions, considering the coefficients of e k 
as functions of 9. Introducing 



u = 4 sin 2 §, 



u = A sin 2 9, A% = ai + a 2 , 
we get the following relations at orders e° and e 1 : 



\ 11 

A 2 = a 1 + a 2 



00 1 ni / 
cotfE^T-f^i-^: 



00 1 {P 

i „^ 1 « j 

■K COt & > -T— T- — 



(Si - S x ) . 



(B.12) 
(B.13) 
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At order e 2 , comparison of the coefficients of A 2 and A\ yields 



cot |E 



1 U> 



1 W 



/ /0 , -s 2 = i <-ot v , , 



S- 2 



(B.14) 



cot |E 



1 u j 



& (?) i 



\ (S 2 + S"fj — — gjL?2 + f'S'l 



icot^E 



1 



5 (S 2 + Sf) - - ±S 2 + iSf] . (B.15) 



All these equations are satisfied by analytic expressions for these sums in terms of 9 from 
Section 2 of this paper. Remember that one needs to substitute 9 — > 29 as an argument on 
the r.h.s. 

At order e 3 , we have two independent structures, A\ and A\A 2 . It is more convenient, 
however, to compare the coefficients of A\ and ^(A^ — A 2 ). The first equation coming from 
the coefficients of A\ yields 



-t § ± 4yT (¥0 + |Ci + |C 2 ) =iootef;iy (^C + | Cl + |C 2 ) , (B.16) 
i =1 \ i I J J =1 { j ) J 



where Cj are the combinations of the harmonic sums defined in Eqs. (|2.5jl . This equation is 
also satisfied, if we use analytic expressions for these sums in terms of 9 given in Section 2. 
The second equation at order e 3 comes from the coefficient of Ai(A^ — A 2 ), 



cotfE7^YTH^ + ^ + ^) = l cot ^ 1 



If we introduce two functions 

tf 2 (0) 



2/ 



00 1 ?/■?' 

cotIE7i7Y-(^ 2 -53) 

cot I £7^7(53 -2S& 
i =1 \ J ) J 



we obtain an interesting relation between them, 

3*i(0) + # 2 (0) = J*i(20) + i* 2 (20) . 



I^i^ + ^S 2 Si + ^53) • 
(B.17) 

(B.18) 
(B.19) 

(B.20) 



This relation is satisfied by the explicit results for the sums involved given in Eqs. (|2.67|) . 



C Further results for the inverse binomial sums 

For completeness, in this Appendix we collect some results for the multiple inverse binomial 
sums of lower weights. They can be obtained by applying the operator u(d/du) to our 
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results presented in sections 2 and 3. These sums occur in lower terms of the e-expansion of 
hypergeometric functions given in Eq. ()2.3())l . 

First of all, we list explicit results for some particular cases of the general formulae (|1.2|) 
and ()2.33j) , in terms of the angular variable (jl.3|) : 



h (?) i 

™ 1 V? 

h (?) ? 
h (?) f 

00 1 u j 

£(!F S2 



tan | , 



ln2 



2 ' 



-2Ls ( l ] (9) + Al e [Cl 3 (9) + 0C1 2 (9) - C 3 ] + 9 2 l 



2i2 



1 /)4 



(C.l) 
(C.2) 
(C.3) 

(C.4) 



Then, let us present a few more complicated examples of the results in terms of 9, 



x 1 Xi? 

i=i 



E T^^^i = tan f i C1 2 (tt - 0) \28L 2 e - 8L e l e + \9 2 ] + Cl 2 (20) [2L 2 + \9 2 



+4L,Ls 3 (0) - 2L,Ls 3 (20) + 4Lsc 2i3 (0) - |0 2 C1 2 (9) ~ Cl 4 (0) 
+ [Ls 3 (tt - 0) - Ls 3 (tt)] [Al e - 28L e ] - 80 [Cl 3 (tt - 0) - Cl 3 (tt)] 
+8 [Ls 4 (tt - 0) - Ls 4 (tt)] - 16 [Cl 4 (tt - 0) - Cl 4 (tt)] 

+\6\ - §0% - 86L 3 e + 13( 3 9 + A9leL 2 e \ 



E 



1 v? 



2/ 



3=1 V j ' J 



5i + ^2 



tanf <^ Cl 2 (vr - 



20 2 + 2ZS - 20L e l e + 32L; 



(C.5) 
|Ls 4 (26) 



+C\ 2 (29) 

+ [Ls 3 (tt - t 
22 

+ y [Ls 4 (TT 



1q2 



9 Z - 2L e k + ALt 



+ Ls 3 (26) [l e - AL e ] - 2Ls 3 (9) [l g - 5L 6 



Ls 3 (tt)] [10le - 32L e ] - 69 [Cl 3 (tt - 9) - Cl 3 (tt)] - Cl 4 (9) 



Ls 4 (tt)] + lOLsc 2 , 3 (0) - 12C1 4 (tt - 9) + §Ls 4 (0) 



-\9 2 C\ 2 (9) + 1OC 3 + \9\ - \9 3 L e - 89L 3 e + 89L 2 9 l e - 29L e l 2 e 



(C.6) 



1 v? 



S? + 35^2 + 2S 3 = tan \\ 9C1 2 (tt 



1 2 + (2L e -/ e 



18Lsc 2 , 3 (0)-|Cl 4 (0)-f0 2 Cl 2 (9) 



+ |C1 2 (20) [|0 2 + (2L e - l e ) 2 

-\8(2L e - k) [Ls 3 (tt - 9) - Ls 3 (tt) - ±Ls 3 (9) + ±Ls 3 (29)\ + f Cs0 
+6 [Ls 4 (tt-0)-Ls 4 (tt)]-60 [Cl 3 (tt-0)-C1 3 (tt)]-Ls 4 (20) + 3Ls 4 (0) 



-12CL (tt- 



±0 3 Z e - \9 3 L e + 80Z 3 , - 60Z^L e + 129l s L 2 9 - 89L 3 - 79l 3 e 



(C.7) 



37 



Using the analytic continuation procedure described in Section 3, these results can be 
rewritten in terms of the conformal variable y defined in Eq. (13. 2 j) . Below we list these 
analytically-continued results for Eqs. (jC.lj) - (jC.4|) . as well as for those sums from Section 2 
whose analytical continuations were not presented in Section 3: 



t 1 ' 



i=i 



J 



2j 

™ 1 vP 

i 1 ul 



3=1 



J' 



23 



3=1 



3 



00 1 v j 

3=1 



°° 1 U j - 
3=1 



(1)> 



si 



00 1 vP 

3=1 



(1)1 



^ 1 ui 
^ 1 

S(f)? 2 



i-y 

l + y 



In y, 



2Li 3 (y) - 2 In yLi 2 (y) - In 2 y ln(l - y) + ± In 3 y - 2( 3 



(C.8) 
(C.9) 
(CIO) 



4S 2 , 2 (y) - 4Li 4 (y) - 4S 1;2 (y) lny + 4Li 3 (y) ln(l - y) + 2Li 3 (y) lny 
-4Li 2 (y) In y ln(l - y) - In 2 y ln 2 (l - y) + § In 3 y ln(l - y) - £ In 4 y 



-41n(l-y)C 3 + 21nyC 3 + 3C4 

1-2/ r 



l + y 

1-2/ 

l + y 



2Li 2 (-y) - 2 In y ln(l + y) + \ In 2 y - & 



Li 2 (y) - 2Li 2 (-y) - 2 In y ln(l + y) + In y ln(l - y) 



(C.ll) 
(C.12) 



+ iln^y-2C 2 

i-y 



(C.13) 



8S 1)2 (-y) - 4Li 3 (-y) + 8Li 2 (-y) ln(l + y) + 41n 2 (l + y) lny 



l + y 

-2 ln(l + y) In 2 y + | In 3 y + 4( 2 m(l + y) - 2( 2 lny - 4( 3 
1-1/ 



(C.14) 



l + y 



10S 1)2 (-y) - S 1)2 y 2 + Li 3 (y) + 2S lj2 (y) - 3Li 3 (-y) 



-2 ln(l - y)Li 2 (-y) + 8 ln(l + y)Li 2 (-y) - 2 ln(l + y)Li 2 (y) 
-21nyln(l - y) ln(l + y) - lnyln 2 (l - y) + 4 lny ln 2 (l + y) 
+ | In 2 y ln(l - y) - § In 2 y ln(l + y) + ^ In 3 y + 6( 2 ln(l + y) 

-C 2 ln(l-y)-|C 2 lny-fC3 , 



1-2/ 



6(1 + 2/) 

i ln4 y > 



ln 3 y 



(C.15) 
(C.16) 

(C.17) 
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00 1 u j 



l + y 

i-y 

l+y 



^ In 4 y + 6Li 4 (y ) + In 2 yU 2 {y) - 2( 3 In y - 4 In yLi 3 (y) - 6(4 



(C.18) 



1 In 3 y ln(l +y) - ± In 4 y + k 2 In 2 2/ + In 2 yU 2 (-y) + In 2 yLi 2 (y) 



00 1 u i 



+C3lny-41nyLi 3 (-y)-41nyLi 3 (y)+C 4 +8Li 4 (-y) +6Li 4 (y)J,(C.19) 
I In 3 y ln( 1 + y ) - | In 3 y ln( 1 - y) - i In 4 y + (2 In 2 y + 4( 3 In y 

(C.20) 



l + y 

+ In 2 yLi 2 (-y) - 4 In yLi 3 (-y) + 8(4 + 8Li 4 (-y) - Li 4 (y) 
For the combination (|2.23|) involving (5*2 + Sf) we get 



£ 7^ + S?) = ^] 12S ll2 (-y) - 2S li2 (y 2 ) + Li :1 ( ,/) + 6S., +//) N< :j 



l + y 



-2Li 3 (-y) - 4 ln(l - y 2 )Li 2 (y) + 2 ln(l - y)Li 2 (-y) + 8 ln(l + y)Li 2 (-y) - 2( 2 In y 
+ lny [ln(l - y) - 2 ln(l + y)f + (± In 2 y - 4( 2 ) [ln(l - y) - 2 ln(l + y)] } , (C.21) 



For a separate term of this sum involving S* 2 , we get 



2Li, 



-co* 



2Li 3 (c^)-^ln 3 y 



(C.22) 



where uj s is defined in Eq. (|3.6j) . 

The results for the multiple inverse binomial sums presented in this paper are summarized 
in Table 1. When two equation numbers are present, the first one refers to the result in terms 
of the angular variable 9, whereas the second one corresponds to its analytical continuation. 
An asterisk means that the corresponding equation holds for general c. The symbol f means 
that the results for the sums involving the combinations Si (§ 2 + S'fj and Sf + 3SiS 2 + 25* 3 
can be extracted from the expressions ()3.10|) and ()3.11|) given for the sums involving C\ and 
C 2 , respectively, using the definitions ()2.4|) . 

D Connection between binomial, harmonic and inverse 
binomial sums 

Using the procedure described in section 2.2, one can also construct the e-expansion of 
hypergeometric functions of the following types: 



p+iFp 3 



§+&i£, . . . , § + &/£, l+ai£, . . . , l + a K e, 2 + dxs, . . . , 2 + d L e 



+he, 



' 2 



+/j-i£, l+ei£, . . . , l + e R e, 2+de, 2 + c K+L - R e 



u 
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Table 1: Equation index for the inverse binomial sums 



p+i^p 



§+&i£, . . . , \+b J E 1 1+aie, . . . , l+a^e, 2+d^, . . . , 2+d L e 



2; 



v |+/i£, . . . , |+/j£, 1+eje, . . . , 1 + 2 + cie, . . . , 2+c^+L-i?-^ 

where P = K + L + J— 1 and w = 4z. By analogy with Eq. ()2.29|) . the e-expansion of these 
functions can be written in the following form: 



-iF F 



il + hsy, {l + a i£ } K , {2 + d i e} L 
{1 + fte} J -\ {1 + e i£ } R , {2 + c t e} K+L - R 



u 



2n 



K+L-R, 
=1 I 



l + c 8 e)Il£ 1 1 (l + 2/ fc e 



-1^ 



w nf =1 (l + di£)Ilt =1 (l + 2b r e) 

{'l + bie} J , {l + ai e} K , {2 + die} 1 
{§ + ^} J , {1 + e,e}* {2 + c^}^-*- 1 

1 nf^-^l + c s e)U{ =1 (l + 2f k e) « 




(D.l) 



2 J 



nf =1 (i + die)n 



(D.2) 



where the function A is defined in the same way as in Eq. ()2.30|) . One should only remember 
that the upper summation limits for the coefficients Bk and Ck are changed, since the 
numbers of the parameters bi and Cj in Eqs. (|D.1J) and (|D.2|) are different. 

The sums appearing in Eq. (jD.l|) are expressible in terms of the multiple binomial sums 
[8], whereas the sums of Eq. (jD.2|) are reduced to the multiple harmonic sums. Using relations 
between hypergeometric function of different arguments, it is possible to express one type 
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of sums in terms of the another one, plus some trivial part. Let us illustrate this on the 
example of 2-F1 function (|2.1jl . Using the standard formula of analytic continuation to the 
argument 1/z, we obtain a combination of two 2 F\ functions. To bring them to the form 
of (ID.lj) , we need to shift some of the parameters using 



1F1 



a, b 

c 



ahz ^ 1 1, a + 1, 6 + 1 
2, c+1 



3-^2 



Finally, using the duplication formula for the argument of the T-function, we arrive at the 
following relation: 



x r(l + 6e) „ 



1 + axe, 1 + a 2 e 
! + be 



r(l + (o 2 - ai)e)r(l + {b - ai)e) 



(_ 4z )a l£ r(l + a 2 e)r(l + 2(6-ai)e) 



f (l + o 1 e)(l+2(o 1 -6)e) 
I + 2s(l + (01 - a 2 )e) 3 2 



1, 2 + ai£, § + (a-i - 6)e 



2, 2 + (ai - a 2 )e 



_L r(l + (ai-a 2 )e)r(l + (6-o 2 )e) 

_4 z )a 2£ r(l + aie)r(l + 2(6 - a 2 )e) 



J (l + a2 £)(l + 2(a 2 -6) £ ) 
\ 2*(1 + (a 2 - ai)e) 3 2 



1,2 + a 2 e, § + (a 2 - 6)e 
2, 2 + (a 2 - ai)e 



(D.3) 



Here, the e-expansion of the hypergeometric function on the l.h.s. can be expressed in terms 
of the multiple inverse binomial sums (see Eq. (I2.29j) ). whereas the functions on the r.h.s. 
yield the multiple binomial sums (see Eq. (jD.lJl ). 

Furthermore, the following three quadratic relations for 2 iq functions (see, e.g., in [31]) 
allow us to connect some multiple binomial sums (appearing in Eq. (jD.ip ) with multiple 
harmonic sums (see Eq. (ID . 2|) ) : 



2 Pi 
2-P1 



3 
2 " 

2 + (a+b)e 
3 

2 + 2be 

2-1 « 



l+6s, | + ae 



l+6e, | + a£ 



(i + x) 2+2a£ 2*1 



(l + X ) 3+2a£ 2*1 



1 



l-2e 



2(l-£l 



2 f x 



2+2ae, l + (o-6)e 
2 + (a+6)£ 

ae, 1 + (a — 6)s 



| + 6e 



( ^ 






•)] 





1 + 2*1 



l,e 
2-e 



(D.4) 
(D.5) 

xj , 

(D.6) 



where 



X 



1 - vT^tx 



u 



4x 



(D.7) 



1 + ~ (l + x) 2 ' 

and all orders of the ^-expansion of the last 2 Fi function on the r.h.s. of Eq. (jD.6|) are known 
through Eq. (2.14) of Ref. [16]. 
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Combining all these relations together with the results for the multiple inverse binomial 
sums presented in this paper, we have reproduced several known results for multiple binomial 
sums, including those obtained in Ref. [8]. For completeness, we list a number of such results, 
including trivial ones, which could be extracted from Ref. [8] (although were not explicitly 
listed there), 



E 

i=i 

oo 

E 

i=i 

oo 

E 

oo 

E 

oo 

E 

i=i 

oo 

E 

i=i 

oo 

E 

i=i 

oo 

E 

i=i 

oo 

E 

i=i 

oo 

E 

i=i 

oo 

E 



E 



.7) •' 

'2 A zj_ 
J J 3 

.3)3 
'2 A z^' 



2j 



2X 

i-x' 

21n(l + X ), 

-2Li 2 (- X )-21n 2 (l + X ), 



,3)3 
'2j 



= 4S 1>2 (- X ) - 2Li 3 (-x) + 4Li 2 (- X ) ln(l + x) + ! ln d (l + x), 



1-X 



[(1-X) ln(l + x) + (l + x) In(l-X)], 



^)y5 1 = Li 2 ( X 2 )+21n 2 (l + x), 



'2j 

'2 A z^' 



2* Si 



[Xln(l + x)-(l + x) In(l-X)], 



,3)3 
'2j 



-S 1= 2Li 2 ( x ) + ln 2 (l + x ), 
4 X T • 



z 3 S 2 = 



2 



Li 2 (- X ) + 21n 2 (l + X ), 



i-x 



i-x 



4 X Li 2 (x) + 2 X Li 2 (-x) - (1-x) ln 2 (l + x) + 2(l + x) ln 2 (l-x) 
6 X Li 2 ( X ) - 2(1 + x) ln(l + X) m(l " x) 



(s 2 - st 



-(1 - x) ln 2 (l + X) + 4(1 + X) ln 2 (l - x) 

f 2 X Li 2 (x) - 2(1 + x) ln(l + X) m(l - x) 
+ X hi 2 (l + x) + 2(l + x)hi 2 (l-x) . 



i-x 



(D.8) 



Furthermore, explicit results for the multiple binomial sums with 1/j 4 , Si/j 2 , Si/j 2 , 
S 2 /j, S 2 /j, SiSjj, and (S 2 - S 2 )/j are presented in Eqs. (A.4)-(A.10) of Ref. [8] (where 
the same notation x as here was used, with x corresponding to our z). We confirm all those 
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results. Moreover, using our approach we have obtained a number of new results for the 
multiple binomial sums, 



jt( 23 ,) z>S 3 = -|ln 3 (l + x )-41n(l + X )Li 2 (- X )- I ^ [2S 1>2 (-*) + X Li 3 (-*)] , 



J 



J Z 

" h,:1(l f llu(:l + \) Li 2("\)+ 1^(1+ v)Li 2 (\) 

j=l V J J 

4 



+ 



i-x 



X Li 3 ( X ) + Si, 2 (x 2 ) - 2Si, 2 (x) + (1 + x) ln(l - x)Li 2 (-*) 



E 



2j 



z^S 2 S 1 = I ln 3 (l + x) + 4 ln(l + X )Li 2 ( X ) 

+ [2 X Li 3 ( X ) - 6S 1)2 (-x) + 2S 1>2 (x 2 ) - 4S 1)2 (x) - xLi 3 (~x) 

, 2(1 + X) 



i-x 



[2m(l-x)-ln(l + x)]Li 2 (-x) , 



E 

j'=i 



'2/ 



i S 1 (S 2 - Sf 



2(1 +X) 

i-x 



= 7^— [16Sl - 2 (x) " 2xLi3 (x) " 6xln(1 + x)Li2 

^ X 

41n 3 (l-x)-41n 2 (l-x) ln(l + x) 



+ ln(l - x) ln 2 (l + X) + 8 ln(l - x)Li 2 (x)] + | ln 3 (l + x) , 



2S 1;2 (x 2 ) - 6Si, 2 (-x) + 16Si, 2 (x) - xLi 3 (-*) 



2(1 + X) 

i-x 



i-x 

-4 X Li 3 (x) + 2(1 - 3x) ln(l + x)Li 2 (x)] 
4 ln 3 (l - x) - 2 ln 2 (l - X ) ln(l + x) + 2 ln(l - x)Li 2 (-*) 
+ 10 ln(l - x)Li 2 (x) - ln(l + x)Li 2 (- X )} ~ § ln 3 (l + x) , 



U1> 



(35i5 2 - 2S 3 - S 3 ) = - 
2(1 + X) 



-x 



12Si )2 (x) -6x ln(l + x)Li 2 (x) "X ln 3 (l + x) 



i-x 



6 ln(l - x)Li 2 (x) - 6 ln 2 (l - X ) ln(l + x) 
+31n(l-x)ln 2 (l + x)+41n 3 (l-x)" , 



'2j 



£ P ) zjS i = -I ln3 (! + X) - 41n(l + x)Li 2 (-*) - 121n(l + x)Li 2 (x) 

3=1 \ J J 

' [3S 1)2 (x 2 ) - 4S 1)2 (-x) + 6S 1)2 (x) - 2 X Li 3 (-*) - S X U 3 ( X ) 



i-x 

4(1 + X) 

i-x 



21n 3 (l-x) + 31n(l-x)Li 2 (-x)+61n(l-x)Li 2 (x) , (D.9) 
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where x is defined in Eq. (|D.7J) and, as before, u = Az. These results were recently used to 
construct the e-expansion of two-loop sunset-type diagrams of special type (for details, see 
section 3.1 in [62]). 

The results for harmonic sums of the type Y%Li z^S ai . . . S^S^ . . . could be deduced 
from the following relation: 

4(a-6)e 



2^1 



1, f+ae 
l + be 



r(2 + 26e)r(l + oe)r(l + (a - b)s) 
r(l + be)Y(2 + 2ae) 



(1 _ z y+(a-b)e z l/2+be 



l + 2be 



2^1 



1, '^ + ae 
2 + (a-b)e 



1 - z 



(D.10) 



2(1 + (a - b)e) 

which corresponds to an analytic continuation of the 2 -Pi function from the argument z to 
(1 — z). In particular, we get 







1 - ^ 
1 



i ln 2 (l - *) + \ ln 2 (l - Vi) + \ ln 2 (l + v^) + §Li 2 



2v^ 
1 



Li 2 



1 + x/i 



- Li 2 



4^ 



ln 2 (l - V5) - ln 2 (l + y/z) 



In 2 
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In 



Introducing a new variable £ such that 



s 



i + V^ (i + 2 

we can present the above result in a more compact form, 



£SiS: 



1^ 



3=1 



2(1-*) 



Li 2 (z) + ln 2 (l-z) 



4£ 



Li 2 (-0+ln£ln2 + |C2 



+ 



(D.ll) 



(D.12) 



^In^ 



;d.i3) 



where we have taken into account that 



Li 2 



1 + x/i 



Lis 



1- v^" 



2Li 2 (-0+m£ln(l + + C2- 



Integrating the representation (|D.13j) . it is easy to get the following result: 



3=1 J 

-2Li 3 (-0 + Li 2 (-0 In £ - |C2 In £ + £ In 3 £ . (D. 14) 

In Appendix E of Ref. [14] the results (up to weight 4) for the harmonic sums of the 
type Y%=i Sai ■ ■ ■ S ak zi / j a have been presented. It is easy to extend those results to the case 
E£li S ai ... S ak zi/j a , using the following property [32]: 



-Si, 2 (z)-i]n(l-z)Li 2 (*)-§ln 3 (l 



C 3 -|ln 2 £ln(l-^) 



j= i J 



then £ f(2j)- = 2 a - 1 [f{^z) + F(-y/z) 
3=1 J 



;d.is) 
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As a consequence, we get 



E^ = (v^) - ^7i) Li « (-v^) > 



■>a— 1 



Sa-1,2 (v 7 ^) + Sa-1,2 (~V^) 



(D.16) 



As a further illustration, we also present results for other sums up to weight 3, 



E^v 

00 7 3 

J'=l ^ 

00 

E^v 
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2(1 + y/z) 



Li 3 (v^) - ln(l - v^)Li 2 (y/Z) - Si )2 (y/Z) 



Lis (-y/Z) - ln(l + v^)Li 2 (-Vi) - S lj2 (-^i)] . (D.17) 



To investigate relations between multiple harmonic sums the technique proposed in Ref. [63] 
can be useful. 

As a illustration, we present the hihger order e-expansion of the some of the hypergeo- 
metric functions: 
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1 + ai£, 1 + a 2 e 
2 + ce 
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-m(l -z) -s< 
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ln 2 (l -z) + cLi 2 (z) 
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ln 3 (l - 2) 



-£ 3 C 
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+ (c - ai)(c - a 2 )(c — ai — a 2 ) 



ln(l-z)S 1)2 (z) + ^ln 2 (l-z)Li 2 (z) 
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U ' " J a2) -ln\l-z) + c{c-a 1 -a 2 ) 2 S 13 {z) + c 3 U 4 {z)\ + 0{e 4 ) I . (D.li 



24 



E The 0(aa s ) corrections to the polarization function 
of neutral gauge bosons in arbitrary dimension 





Figure 2: Two-loop contributions to the off-shell polarization function of a neutral gauge 
boson. Bold and thin lines correspond to the massive quark propagator and the massless 
boson (gluon or photon) propagator, respectively. 



Here we present an example of a physically relevant calculation 11 that can be expressed 
in terms of the master integrals Jon studied in Section 4.3. Let us consider the two-loop 
propagator-type diagrams shown in Fig. |2j All of these 0(aa s ) contributions to the polar- 
ization function of the gauge bosons involve a quark loop with a gluon exchange. It was 
analytically calculated in [56], up to the finite term of the e-expansion. Here we present 
the bare two-loop results in n-dimensional space-time (see also in Ref. [8]). In contrast to 
the calculations performed in [56], here we use Tarasov's recurrence relations [38] for the 
reduction of the original integrals to the set of master integrals. In this Appendix we use 
the Euclidean notation [57], P 2 <-> — p 2 , so that the on-shell limit would read P 2 — > — m 2 . 

Let us decompose the polarization tensor into the transverse Hr(P 2 ) and longitudinal 
n L (P 2 ) parts, 



n^(P 2 



P 2 



Ht(p 2 ) + 



PnP 



H L (P 2 



Then, the two-loop corrections can be written as 
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11 Another non-trivial example where these master integrals appear is given in Ref. [55]. 
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+8Am 2 tn 3 - 72Am 2 tn 2 + 232Am 2 tn + ASm 2 tn 3 - AASmHn 2 + 152Sm 2 tn 

-32At 2 - 2A0Am 2 t - 32St 2 - 160,Sm 2 t] | , 

n[ 2) (P 2 ) ^^^{-, on( i, 1,1)1^^ 

-Jon(l,l,2 
+ [A (m 



(E.l) 



16m 



4m 2 n 2 - 28m^n + 48m 2 + tn z - 5tn + 8t 



- [B (m,m,t)} 



t(n-A)(n- 3) 
]2 8(n-2) 
J t (Am 2 + t) (n - 4) 
2 8m 2 



— A (m)B (m, m, t) 



(Am 2 + t)(n- 4) 
fi <n - 2) (n 2 



tn 2 — 3tn + Am 2 n — 16m 2 



-Atn + At + Am 2 n 2 - 2Am 2 n + tn 2 + 40m 2 



5n + 8) 



-2m 2 n + 10m 2 + t 



(E.2) 



(4m 2 + t) (n-A)(n-3) 
where the vertices Vj (j = 1, 2) are defined as 12 

v 3 = ' l 91»(vj + %75) , 

and we have introduced the following notations: 

A = aia 2 , S = V\V2 , t = P 2 , 

N c is a color factor (equal to 3 for quark and 1 for lepton), Cp is the Casimir operator of the 
fundamental representation of the Lie algebra (equal to | for SU(3) and 1 for QED), and m 



12 



Explicit values of coefficients Vj and dj for the Standard Model can be extracted from [57]. 
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is the mass of the loop fermion. Finally, the occurring integrals are defined as 

7T-" f r d n Ki d n K 2 



Jon 0,^1,^2) 
B (m a ,m b ,t) 
A (m) 



T 2 3 



7T 



-n/2 



r (3 - f 



[(K 2 - P) 2 ] a [K 2 + m 2 p [(Kx - K 2 ) 2 + m 2 P ' 
d n /\~ 



[(K - P) 2 + m 2 ] [K 2 + m 2 ] ' 



7T 



-n/2 



d n /\~ 



4(m 



2\n-2 



r (3 - |) J K 2 + m 2 (n- 2)(n - 4) 



(E.3) 



We note that the integral Jon is defined in the same way as in Section 4.3, one should only 
remember to substitute p 2 — > —t. 

In particular, for the zero momentum transfer, P 2 = t = 0, we get 



n[ 2) (o) 
ni 2) (o) 



-aM^N c C f [Mm)f 4 (n-2)(n»-5n + 7) 



'(4vr)™/2 



, 9 2 9 2 s Mr U( „ 2 A(n - 2 )(n 3 - 6n 2 + 13n - 11) 
-A^— —N C C F [A (m)\ 



(4tt)"/ 2 



(n — l)m 2 



To obtain the finite terms of the e-expansion of the results given in Eqs. (jE.lj) and 
(IE.2|) . the integral Jon(l, 1, 2) should be expanded up to the e-part. Furthermore, using the 
approach of Ref . [38] , we obtain the following relations between the integrals investigated in 
Section 4.3 and the master integrals Jon(l, 1, 1) and Jon(l, 1, 2): 



(3n-8)(3n-10)(n-3) 2 



+Jon(l,2,2) 



n - 4 
im 2 (t + m 2 ) 
n-A 



[A (m)] 2 (n - 3)(n - 2) : 



(7n - 24) 



4m 4 2m 2 (n - 4) 
- 2t 2 (n - 3) + 2tm 2 (7n - 17) + 8m 4 (2n - 5) 



Jon(l, 1, 2 



[J ou (l, 2, 2) + 2J n(2, 1, 2)] (n - 3)(t + 4m 2 )(t - 2m 2 ) , 
(3n-8)(3n- 10)(n-3) 2 



(E.4) 



+Jon(l,2,2) 



n-A 
A(t + 2m 2 
n-A 



n 



(n-3)(2n-7)(3n-8)(n-2) 2 
~ [Mm)] Am\n-A) 

(ft + 12m 2 ) + n (f * + 62m 2 ) - 55t - 86m 2 



+| [Jou(l, 2, 2) + 2 Jon (2, 1, 2)] (n - 3)(3n - 8)(t + 4m 2 



(E.5) 



From these relations we see that, in order to obtain the e J terms of Jon(l, 1, 1) and Jon(l, 1, 2), 
the integral Jon(l, 2, 2) should be expanded up to e> , whereas the combination [J n(l, 2, 2) + 
2Jon(2, 1, 2)] up to £ J " -1 only. The expansion of the integrals near the threshold can be per- 
formed by using the numerical algorithm described in Ref. [58]. 

The higher-order terms of the ^-expansion of the one-loop integral B (rrii,m 2 ,t) can be 
extracted from Refs. [16,30]. Here we present coefficients up to the order e 2 for the particular 
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case mi = m 2 = m: 



B (m, m, t) 



m 



2\-e 



+ 



1+y 



;i-2e)Ve 1-y 



my + e \ ln 2 ?/-21nyln(l + y)-2Li 2 (-y)-( 2 



4Si, 2 (-y) - 2Li 3 (-y) + 4 ln(l + y)Li 2 (-y) - In 2 y ln(l + y) 



+2 lny ln 2 (l + y) + ± In 3 y - ( 2 In y + 2( 2 ln(l + y) - 2( 3 + 0(e 



(E.6) 



where the variable y is defined in Section 3. 

The bare two-loop amplitudes (jE.l|) and ()E.2|) contain sub divergencies which should be 
canceled by proper counterterms, 

nc T (P 2 ) = W^n«(p 2 ), 

where is the one-loop amplitude, while 5m 2 is the one-loop mass counterterm defined 
in a particular renormalization scheme. The derivatives of the bare one-loop amplitudes in 
n dimensions read 

,2 



9 :n«fp 2 ) 



dm 2 T 







dm 2 L 



Ii) l) (P 2 ) 



A 



(4^)n/: 



:iV c 4< 



St(3-n) 
(Am 2 +t) 
t(n - 3) 



(4m 2 + 1) 



^S+(2-n)A 
+ l-n 



x 25(71-2) , , , 
B (m, m, t) + — — -— A (m) 
{Am A +t) 



2(2 -n) , . J 



As an example of application of these formulae, let us consider the transversal part of 
the 7 — Z propagator (A = 0) in the MS'-scheme. Up to the finite in e part, the result for 
the subtracted quantity n sub (P 2 ) = lT 2 )(P 2 ) + n CT is 



-rrsub / p2\ n 9 9; 



(4tt)- { e 
<l-4y + y 2 ) 



1-4- 



-fm 2 



v0--v¥ 



i-y)\ 

[ln(l - y) + 2 ln(l + y)] lny [(1 + y 2 ) lny - 2(1 - y 2 ) 



8 2 (2 + 7y-22y 2 + 6y 3 ) 1 2 2 (1 - 6y - 46y 2 - 6y 3 + y 4 ) 

+|m - — In y — 4m lny- 



+4t In ■ 



m 



(i-y) 2 

12y 



y(i-y l 



'i-yf 



96m 2 - — — — — In y In — - 



i-y 2 



[Li 2 (y) + 2Li 2 (-y)] [l - y 2 - 2 (l + y 2 ) lny 



,32^2 (1-^ + y 2 
3 Z/(l-Z/) 2 

+32m 2 (1 ~ 4 ^f ) (l + y 2 ) [Li 3 (y) + 2Li 3 (-y)] + 0(e) 
where we have taken into account the one-loop massive counterterm 

5m 2 = _ 2 jLc F -m 2 . 



(E.7) 



(4tt) 
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The small- momentum expansion of Eq. (|E.7j) is 

^tUp 2 )\ p ^ = S ^ 2 N c C F {lu + fu + + In ^ (An + §u 2 ) + 0(u* 



p 2 ^o (Air) 2 c F \ e 3 405 

where u = —P 2 /m 2 . Let us remind that in a theory with the spontaneous symmetry 
breaking, like the Standard Model, the inclusion of tadpoles [59] is also required for the 
renormalization group invariance of the massive parameters [8,58,60]. The proper bare 
two- loop tadpole contribution is given in Section 4.3 of Ref. [8]. 
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